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ELASTIC  STRESS  ANALYSIS  OF  AN  ADHESIVE  LAP  JOINT 
SUBJECTED  TO  TENSION,  SHEAR  FORCE  AND  BENDING  MOMENTS 


Murlidhar  H.  Pahoja,  Ph.D. 

Department  of  Theoretical  and  Applied  Mechanics 
University  of  Illinois  at  Urbana- Champaign,  1972 

A  stress  analysis  of  the  lap  joint  is  presented  treating  the  problem  as  one  of 
plane  strain.  The  joint  is  subjected  to  a  general  loading,  consisting  of  tension,  shear 
force  and  bending  moments.  The  variation  in  the  material  properties  and  thickness 
of  the  two  adherends  is  considered.  The  displacement  field  in  the  adhesive  layer  is 
expressed  in  series  form  and  the  compatibility  condition  at  the  interface  is  used  to 
express  the  displacement  field  in  the  adherends.  The  potential  energy  of  the  joint  is 
calculated  and  minimized  to  obtain  linear,  ordinary  differential  equations  and  boundary 
conditions.  The  differential  equations  are  solved  on  the  computer. 

Photoelasticity  is  used  to  confirm  the  theory.  Two  specimens  of  lap  joint  using 
a  1/4  in.  layer  of  a  photoelastic  plastic  simulating  the  adhesive,  are  tested  photoelas- 
tically.  The  agreement  between  the  theoretical  and  the  experimental  results  is  found 
to  be  good. 

Design  recommendations  for  the  lap  joint  are  made  based  on  the  results  of  this 


investigation. 
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1.  INTRODUCTION 

In  the  last  two  decades,  adhesive  bonding  has  become  more  and  more  common 
in  engineering  structures.  The  advantages  of  the  adhesive  joint  over  the  conventional 
mechanical  fasteners  are,  savings  in  weight  and  cost  and  elimination  of  holes  which 
can  cause  excessive  stress  concentration. 

The  commonest  type  of  joint  in  use  today,  is  the  lap  joint  where  two  adherends 
overlap  and  are  fixed  together  by  a  layer  of  adhesive  between  them  (Fig.  la).  The 
main  advantage  of  this  joint  is  its  simplicity  in  design.  Further,  where  the  surfaces 
need  chemical  preparation,  the  adherends  are  easy  to  handle.  It  is  also  a  simple 
type  of  joint  to  cure.  The  chief  disadvantage  of  this  type  of  joint  is  that  when  loaded 
in  tension,  varying  shearing  and  tearing  stresses  are  produced  in  the  adhesive  as  well 
as  the  adherends,  and  very  high  peak  stresses  are  produced  at  or  very  near  the  ends 
of  the  overlap. 

Failure  in  a  joint  may  occur  either  in  the  adherend  or  in  the  adhesive  or  at 
the  interface.  A  failure  at  the  interface  is  termed  as  'adhesive, '  since  it  involves 
the  failure  of -the  bond  between  the  adhesive  and  the  adherend.  When  the  fracture 
occurs  either  inside  the  adhesive  or  inside  the  adherend,  the  failure  is  then  termed 
as  'cohesive. '  In  a  cohesive  failure,  the  material  at  the  two  fracture  surfaces  is  the 
same.  The  interfacial  adhesion  of  modem  adhesives  is  high  so  that  failure  mostly 
occurs  either  in  the  adhesive  as  in  the  case  of  metal  to  metal  joints,  or  in  the  ad¬ 
herends  as  happens  with  wood  joints. 

The  stress  in  the  adhesive  layer  is  dependent  on  the  geometrical  proportions 
and  the  elastic  constants  of  the  materials  of  the  joint,  and  on  the  nature  of  the  forces 
acting  on  the  joint.  When  loaded  in  tension,  non-uniform  shear  and  tearing  stresses 
are  developed  in  the  adhesive  layer.  The  shear  stress  is  parallel  to  the  interface 


2 


while  the  tearing  stresses  are  normal  tc  the  in  e _ face.  It  has  been  predicted  theoreti¬ 
cally  and  confirmed  experimentally  that  the  pe.  ..  shear. and  tearing  stresses  occur  at 
‘hr*  two  ends  of  the  overlap.  This  concentration  of  stresses  at  the  ends  of  the  overlap 

i 

is  due  to  two  main  causes,  namely  to  the  differen.ial  strain  of  the  adherends  and  their 
bending. 

I 

Fig.  la  shows  a  lap  joint  loaded  in  tension.  Each  adherend  bears  the  full  load 
P  just  before  the  joint  and  transmit  it  gradually  to  the  other  through  the  adhesive. 

Thus  the  stress  in  adherend  (I)  will  be  the  highest  at  A  and  gradually  dimmish  towards 
B  wh  ?re  it  will  be  zero.  On  the  contrary  the  stress  in  adherend  (II)  will  be  the  highest 
at  B  and  diminish  to  zero  at  A.  The  adherends  (I)  and  (II)  will  develop  tensile  strains 
proportional  to  the  stiesses  resulting  in  the  deformation  as  shown  in  Fig.  lb.  The  dis¬ 
placement  is  the  highest  at  the  two  ends  oi  the  overlap  and  therefore  the  highest  shear 

i  > 

stress  in  the  adhesive  can  be  expected  to  occur  at  these  two  points. 

As  can  be  seen  in  Fig.  la,  the  two  equal  and  opposite  tensile  forces  acting  on 
the  joint  are  not  colinear,  and  the  joint  therefore,  is  not  in  equilibrium.  In  order  to 
overcome  this  difficulty,  it  is  generally  assumed  that  the  line  of  action  of  the  two 
forces  is  oblique  arid  passes  through  the  midpoint  of  the  joint.  Tins  is  equivalent  to 
the  addition  of  a  shear  force  V  (Fig.  2)  so  that  the  moment  due  to  P  is  balanced  by 
the  moment  due  to  V.  This  configuration  of  forces  produces  bending  along  the  joint. 

The  stresses  developing  in  this  way  are  normal  to  the  interface,  and  are  called  "tearing" 
stresses;  they  are  highest  at  the  areas  adjoining  the  ends  of  the  overlap. 

The  earliest  theoretical  analysis  of  lap  joints  is  that  done  by  Volkersen  (1), 
in  which  the  stresses  arising  from  the  differential  straimng  are  considered  but  the 
tearing  stresses  are  ignored.  The  next  analysis  is  that  of  Goland  and  Reissner  (2) 
who  took  the  bending  deformation  of  the  adherends  into  account  as  well  as  the  tearing 
stresses  in  the  adhesive.  Subsequent  modifications  to  the  theories  of  Voikexsen  and 
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Goiand  and  Reissner  have  been  made  by  Plantema  (3)  and  Kelsey  and  Benson  (4),  These 
theories  are  reviewed  and  discussed  in  detail  in  Chapter  IL 

The  results  from  these  theories  indicate  a  non-uniform  distribution  of  shear 
stress  with  peaks  at  the  endo  ot  the  overlap.  While  in  a  qualitative  sense  this  distri¬ 
bution  is  m  agreement  with  what  has  been  observed  experimentally,  the  magnitude  of 
the  shear  stress  in  the  vicinity  of  the  ends  of  the  overlap  is  found  to  be  in  error.  This 
is  a  serious  shortcoming  since  the  highest  stresses  occur  in  the  same  area.  This 
inaccuracy  is  mainly  due  to  the  assumption  made  in  these  analyses  that  the  stresses 
do  not  vary' m  the  direction  of  the  thickness  of  the  adhesive.  Irrespective  of  how  small 
the  thickness  of  the  adhesive  layer  may  be,  it  has  to  be  taken  under  consideration  for 
the  stress  distribution  in  the  areas  distant  by  the  same  order  of  magnitude  from  the 
edge.  Another  important  limitation  of  the  above  theories  is  that  the  joint  is  considered 
to  be  loaded  in  tension  only.  In  a  general  case,  the  lap  joint  could  be  subject  to  shear 
forces  and  bending  moments  as  shown  in  Fig.  3,  the  only  restriction  being  that  the  ex¬ 
ternal  forces  and  bending  moments  be  in  equilibrium. 

The  objective  of  the  present  w'ork  is  to  analyze  the  lap  joint  when  subjected 
to  the  general  loading  (Fig,  3)  and  to  predict  the  stresses  more  accurately  by  allowing 
them  to  vary  through  the  thickness  of  the  adhesive  layer.  The  analysis  does  not  put 
any  restriction  on  the  thickness  of  the  adhesive  and  is  therefore,  valid  for  thin  as  well 
as  thick  adhesives.  The  thickness  of  the  adhesive  however,  is  required  to  be  much 
smaller  than  the  length  of  the  overlap  which  of  course,  is  realistic.  One  major  limi¬ 
tation  of  the  previous  theories  as  well  as  the  present  analysis  is  the  assumption  that 
the  adhesive  is  isotropic  and  linearly  elastic.  The  behavior  of  most  adhesives  is  visco¬ 
elastic  depending  on  time,  temperature  and  previous  history.  Even  when  an  adhesive 
can  be  considered  to  be  elastic,  it  is  more  likely  that  the  stress-strain  curve  is  a  non 
linear  one.  Furthermore,  as  the  adhesive  sets  in  between  the  adherends,  it  may  do 
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so  in  an  uneven  fashion  so  that  it  would  be  strictly  inadmissible  to  consider  the  adhesive 
layer  as  homogeneous  and  isotropic.  The  nature  of  such  variations  is  however,  uncer¬ 
tain  and  any  attempt  to  incorporate  the  viscoelasticity  or  the  non -homogeneity  of  the 
adhesive  into  the  analysis  would  render  the  analysis  extremely  complicated.  Further, 
any  plastic  flow  of  the  adhesive  would  have  a  relieving  effect  on  the  peak  stresses  in  the 
adhesive  and  an  elastic  analysis  would  therefore,  provide  a  built-in  factor  of  safety. 

It  is  felt  that  once  a  reasonable  analysis  based  on  linear  elasticity  of  the  adhesive  is 
available  for  comparison,  it  would  then  be  possible  to  use  such  tools  as  finite  element 
analysis  to  incorporate  the  non-linear  elastic  or  plastic  behavior  of  the  adhesive. 

As  regards  the  adherends,  they  are  assumed  to  behave  as  beams  or  in  other 
words  the  lateral  stress  in  the  adherends  is  considered  negligible  and  their  mechani¬ 
cal  behavior  is  assumed  to  oe  isotropic  and  linearly  elastic.  While  in  the  case  of 
metal  adherends  this  assumption  holds  good,  it  is  not  strictly  valid  for  wood  and  plastic 
adherends.  Wood,  although  linearly  elastic,  is  anisotropic  and  its  mechanical  prop¬ 
erties  vary  with  the  direction  relative  to  the  grain  and  the  way  of  loading.  Yet  for  the 
relative  comparison  between  similar  joints  the  theoretical  results  are  expected  to  be 
sufficient,  since  the  nature  and  directions  of  the  developing  stresses  do  not  vary  con¬ 
siderably  from  joint  to  joint.  In  the  case  of  adherends  which  exhibit  plastic  behavior, 
the  theoretical  results  may  not  be  applicable  since  the  yield  in  the  adherends  will 
greatly  increase  the  differential  strain  and  consequently  the  shear  stress  in  the  ad¬ 
hesive  layer. 

The  analysis  is  based  on  the  theory  of  elasticity  treating  the  joint  as  a  plane 
strain  problem  and  therefore,  it  goes  without  saying  that  the  results  are  valid  only 
for  small  deformations,  and  for  joints  with  sufficiently  large  widths. 

The  experimental  part  of  the  present  investigation  is  designed  to  check  and  con¬ 
firm  the  theoretical  results,  both  along  the  length  of  the  overlap  as  well  as  through  the 
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thickness  of  the  adhesive.  The  stresses  in  the  adhesive  layer  are  determined  photo- 
elastically  since  the  photoelastic  methods  provide  a  simple  and  effective  means  for 
handling  plane  problems  such  as  the  one  tinder  consideration. 
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2.  REVIEW  OF  PREVIOUS  THEORETICAL  AND  EXPERIMENTAL  WORK 


2. 1  Introduction 

Several  papers  have  been  published  analyzing  stresses  in  various  types  of 
adhesive  joints.  In  the  case  of  double  lap  joints  and  scarf  joints  the  analysis  is  rela¬ 
tively  simple  because  of  the  inherent  symmetry  of  the  joints.  The  single  lap  joint 
presents  a  special  case  where  the  eccentricity  in  loading  produces  a  bending  moment 
which  has  to  be  accounted  for  in  the  stress  analysis.  Theories  specifically  dealing 
with  the  single  lap  joint  will  therefore,  be  reviewed  here.  Theories  dealing  with  joints 
such  as  double  lap  joint,  scarf  joint,  butt  joint  or  tubular  joints  will  not  be  included  in 
the  review. 

Experimental  works  have  been  published  describing  tests  to  determine  the 
strength  of  lap  joints.  Such  works  will  be  referred  to  when  useful  in  demonstrating 
a  point  in  relation  to  a  theory.  Of  greater  interest  are  experimental  works  where  the 
main  interest  is  the  distribution  stresses  in  the  joint.  These  works  will  be  reviewed 
in  some  detail. 


2.2  Volkersen's  Theory 

The  simplest  approximate  theory  is  that  attributed  to  Volkersen  (1),  This 

analysis  is  concerned  only  with  the  stresses  arising  from  the  differential  straining 

in  lap  joints  and  does  not  examine  the  tearing  stresses  resulting  from  the  bending  of 

the  adherends.  The  bonded  members  are  assumed  to  be  in  pure  tension  and  the  cion 

gation  of  the  adherends  c  and  the  shear  deformations  of  the  adhesive  7  are  taken 
°  x  xy 

into  account  as  in  shear  lag  problems.  It  is  further  assumed  that  ^  in  the  adhesive 

is  zero  so  that  the  shear  stress  t  is  constant  over  the  thickness  of  the  adhesive 

xy 

layer.  Tire  results  indicate  a  non-uniform  distribution  of  shear  stress  with  peaks  at 
the  ends  of  the  overlap.  Volkersen  compared  the  maximum  shearing  stress  at  the  end 
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of  the  overlap  with  the  mean  stress  and  evaluated  the  stress  concentration  factor  h, 
which  is  given  by. 


n  =  \Fd/W 


W  -  1  +  cosh  \TDW~ 
sinh  n/DW 


where, 

2 

D  =  Ga  /E2t2t 
W  =  (Ejtj  +  E^/E^ 

Here  E^^  >  E2t2>  and  the  maximum  stress  occurs  at  the  interface  with  adherend  (11), 
For  identical  adherends  W  reduces  to  2  and  the  factor  of  stress  concentration 
can  be  written  as, 

n  =  \Td72  coth  n/dTT 

Thus  it  is  seen  that  for  identical  adherends  the  factor  h  of  stress  concentration  is  a 
function  of  a  single  dimensionless  coefficient  D, 

In  the  case  of  a  joint  between  adherends  that  are  not  identical,  the  values  of 
D  and  W  are  different  at  the  two  interfaces.  At  the  interface  with  an  adherend  of 
higher  rigidity,  D  is  smaller  and  W  larger  than  the  corresponding  values  of  D  and 
W  at  the  interface  with  less  rigid  member.  Mylonas  and  deBruyne  (5)  show  in  the 
form  of  a  graph  the  variation  of  Volkersen’s  stress  concentration  factor  h  in  terms 
of  D  for  various  values  of  W.  It  can  be  readily  seen  from  this  graph  that  n  is  smaller 
for  lower  values  of  D  and  higher  values  of  W.  Thus  the  stress  concentration  would  be 
higher  at  the  interface  with  the  less  rigid  member.  The  same  result  is  shown  by 
Greenwood  (6)  by  calculating  the  factors  of  stress  concentration  at  the  two  interfaces 


of  a  joint  between  steel  and  rubber  adherends.  The  conclusion  then  would  be  that  it  is 
disadvantageous  to  have  joints  between  adherends  that  greatly  differ  in  rigidity  since 
this  weald  result  in  a  much  greater  stress  concentration  at  one  interface. 

Volkersen's  theory,  while  it  indicates  in  a  general  way  the  distribution  of  shear 
stress  in  the  adhesive  layer  and  the  importance  of  various  parameters  which  have  an 
effect  on  this  distribution,  has  the  following  shortcomings. 

(1)  Bending  deformation  of  the  adherends  and  the  associated  tearing  stresses  in 
the  adhesive  are  ignored. 

(2)  The  assumption  that  the  stress  does  not  vary  through  the  thickness  of  the 
adhesive  is  shown  to  be  untrue  in  the  proximity  of  the  ends  of  the  overlap,  by  Mylonas  (7) 
and  by  Tuzi  and  Shimada  (8)  in  their  photoelastic  experiments.  This  assumption  results 
in  inaccurate  prediction  of  stresses  in  the  area  where  the  maximum  stress  occurs,  and 
the  factor  of  stress  concentration  obtained  by  this  theory  is  therefore,  unreliable. 

(3)  The  analysis  limits  itself  to  the  case  when  the  members  are  subjected  to 
tensile  loads  only.  Shear  force  and  bending  loads  are  not  considered. 

2.3  Theory  of  Gola nd  and  Reissner  (2) 

This  theory  takes  into  account  the  bending  deformation  of  the  adherends  and 
also  the  tearing  stresses  in  the  adhesive.  The  system  considered  is  that  shown  in 
Fig.  2a.  The  line  of  action  of  the  tensile  force  R  is  assumed  to  pass  through  the  mid¬ 
point  of  the  adhesive  layer.  This  is  equivalent  to  a  system  where  a  small  shear  force 
V  acts  normal  to  P  so  that  the  couple  due  to  P  is  balanced  by  the  couple  due  to  V. 

The  adherends  are  identical  in  terms  of  mechanical  properties  as  well  as  geometry, 

Goland  and  Reissner  recognized  that  the  bending  of  the  adherends  outside  the 
joint  region  has  a  pronounced  effect  upon  the  stress  distribution  in  the  joint  itself.  This 
effect  is  expressed  via  their  bending  moment  factor  k  which  is  obtained  by  treating  the 
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adherends  as  cylindrically  bent  plates.  The  value  of  k  is  given  by: 


=  1  +  2\p2  tanh 


76(1  -  Pj2  >  •  PAEjtp 


The  system  can  now  be  reduced  to  that  shown  in  Fig.  2b  where, 
kPt 

M  --21 
P  =  R 


and 


V  =  kpj3(l  -  U[2)  •  P/(E1t1) 


Goland  and  Reissner  considered  two  cases. 

(1)  In  the  first  approximation,  where  t/t^  «  E/E^,  the  presence  of  the  ad¬ 
hesive  layer  is  ignored  and  the  joint  is  assumed  to  be  a  homogeneous  isotropic  rec¬ 
tangular  slab.  The  calculated  stresses  are  those  of  the  adherends  along  the  glue  line, 
and  are  assumed  to  be  equal  to  those  of  the  adhesive  layer.  The  values  of  the  shear 
stress  t  and  the  tearing  stress  a  are  given  as  converging  series.  The  results  show 
that  while  the  tearing  stress  is  high  at  the  edge  of  the  joint,  the  shear  stress  is  zero. 
The  shear  stress  however,  rises  rapidly  to  a  maximum  in  close  proximity  to  the  edge, 

(2)  In  the  second  case  where  t/tj  »  E/E^  ,  the  adhesive  layer  is  considered 
to  be  flexible  so  that  the  transverse  normal  strain  and  the  shear  strain  in  the  adherends 
may  be  neglected  in  comparison  with  the  corresponding  strains  in  the  adhesive.  In 
this  approximation  the  adhesive  is  assumed  to  behave  like  a  system  of  infinitesimal 
springs  placed  between  two  plates.  The  distributions  of  shear  stress  i  and  the  tearing 
stress  o  are  obtained  as  follows: 


f  •  <Y>  =  T  [4^<1+31:>- 


Ox/t^) 


sinh  (pa  > 


+  3(1  -  k) 


—  ( »  - )  =  — -  (R-o^  •  k/2  +  ^k.  •  cosh  A  •  cos  A)  •  cosh  (-2  — ) 

P  Zl»  Za  Z  I  3. 


.  C0S(~^)  +  (RjX 


k/2  +  ^k^  •  sinhA  •  sinA)  •  sinh  (---—) 


.  , 2\xN 
sin(— ) 


where. 


P  Ejt 


6Et,  i 
A  =  _£__  ( _ L)  4 

2tj  v  Ejt  ’ 


Rj  =  cosh  A  •  sinA  +  sinhA  •  cos  A 


R2  -  sinhA  •  cosA  -  coshA  •  sinA 


A  =  5  ( sinh  2  A  +  sin  2  A  ) 


kj  =  Va/(2Pt1) 


P  =  PAj 
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The  maximum  values  of  t  and  o  are  found  at  the  edge  of  the  joint  and  are  given  by. 


t 

max 

P 


(f)  =  i 
h 


-!%-  (1  +  3k)  coth-||  +  3(1  -  k) 


and 


T 

max 

P 


2 

X  k  #  sinh2X  -  sin2X  . .  cosh2X  +  cos2X 
2  sinh2X  +  sin2X  K1  sinh2X  +  sin2X 


The  first  expression  gives  the  ratio  of  maximum  shear  stress  to  mean  shear 
stress  which  of  course,  is  the  factor  of  stress  concentration.  As  is  seen,  this  factor 
is  a  function  of  /3a/2t^  and  k.  Comparing  with  the  factor  D  of  Volkersen's  theory,  it 
can  be  readily  seen  that  /3a/2t^  =  si  2D.  As  in  Volkensen's  analysis,  the  shear  stress 
distribution  obtained  through  Goland  and  Reis  seer's  theory  is  uniform  irrespective  of 
the  value  of  D  as  long  as  this  is  smaller  than  about  0. 1,  but  the  stress  concentration 
increases  rapidly  for  higher  values  of  D  and  for  D  ^  2  becomes  practically  propor¬ 
tional  to  nTD.  This  is  because  for  large  values  of  D,  coth\?2D  approaches  the  value  1 
and  in  fact  the  factor  ii  of  stress  concentration  can  then  by  written  as, 

\  (1  +  3k)  sf2D  +  3(1  -  k)  / 4 

The  bending  of  the  members  (k  small)  is  seen  to  have  a  reducing  effect  on  the  highest 
stresses.  For  small  values  of  k,  the  shear  stress  distribution  is  of  similar  form  to 
that  obtained  by  Volkensen.  in  the  limiting  case  of  k  =  1  i.  e.  when  the  bending  defor¬ 
mation  of  the  adherends  is  negligible,  the  Goland  and  Reissner  theory  yields  if  -  s[2D 
coth  \r2D,  whereas  the  Volkersen  result  was  ii  =  si  D/2x  coth  \ !  D/2  .  For  D-2  i.e. 
for  long  overlaps,  the  stress  concentration  factor  is  therefore  twice  as  large  as  pre¬ 
dicted  by  Volkersen. 


S5»6^;iBM6i«a^ 


w- 

s 
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The  theory  of  Goland  and  Reissner,  though  an  improvement  over  Volkersen's 
analysis  in  that  it  takes  into  account  the  bending  deformation  of  the  adherends  and  the 
tearing  stresses  in  the  adhesive,  has  certain  limitations  which  can  be  summarized  as 
follows. 

(1)  This  theory  is  valid  only  for  adherends  of  the  same  material  and  of  identical 
length  and  thickness. 

(2)  The  first  case  is  valid  for  t/t^  <  <  E/E^  while  the  second  case  requires 
the  condition  t/t^  >  >  E/E^  . 

(3)  The  joint  edge  loads  are  not  in  equilibrium  except  when  k  =  1. 

(4)  The  stresses  in  the  adhesive  are  not  considered  to  vary  through  its  thickness 
which  as  in  the  case  of  Volkersen's  theory,  results  in  inaccurate  prediction  of  stresses 

in  the  area  adjacent  to  the  overlap  edge. 

(5)  The  theory  does  not  consider  external  shear  force  or  bending  loads. 

2.4  Plantema's  Modification 

Plantema  (3)  has  attempted  to  combine  the  Goland  and  Reissner  theory  with  the 
Volkersen  theory.  Volkersen's  theory  is  employed  to  calculate  the  differential  strain 
of  the  members  and  the  stress  distribution  at  the  edges  of  the  overlap.  The  factor  k 
of  Goland  and  Reissner  theory  is  then  introduced  to  calculate  the  differential  strain 
of  the  members  due  to  bending.  The  total  strain  of  the  members  due  to  both  tension 
and  bending  is  used  for  the  correction  of  Volkersen's  formula.  The  factor  of  stress 
concentration  is  given  by, 

n  =  *VD(1  +  3k)  /2  •  coth  \lD{l+  3k)  /  2 

No  calculation  is  made  of  the  tearing  stresses. 

Mylonas  and  deBruyne  (5)  show  that  for  inflexible  adherends  (k  -  1)  and  long 
joints,  Plantema's  results  approach  those  of  Goland  and  Reissner  whereas  for  flexible 
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adhs rends  (k  «  1)  Plantema's  results  are  closer  to  those  of  Volkersen. 

2.  5  Modification  by  Kelsey  and  Benson  (4,  9) 

Kelsey  and  Benson  employ  the  complimentary  energy  method  to  determine  the 
shear  and  tearing  stresses  in  a  lap  joint.  In  this  method,  the  equilibrium  is  taken  into 
j  account,  but  the  necessary  compatibility  is  only  approximately  satisfied.  Treatment 

*  of  the  lap  joint  as  a  shear  lag  problem  results  in  a  second  order  differential  equation, 

| 

|  but  allowance  for  the  complimentary  energy  due  to  the  tearing  stresses  inflates  the 

|  equation  to  one  of  fourth  order.  The  boundary  condition  that  the  shear  stress  vanishes 

f  at  the  two  ends  of  the  overlap  is  then  applied.  The  expressions  for  shear  and  tearing 

stress  are  obtained  as  follows, 

I 

\  n^a  m.a 

,  (cosh  mjX)  •  (cosh  — j- )  -  (cosh  — ^—)  »  (cosh  m^x) 

i  t  =  5Pm^m2a  •  - 

m.a  m2a  m.a  nua 

(mj  cosh  — 2'  )  •  (sinh  —j—)  -  (m2  sinh  ■  ■  1  )  .  (cosh -tj—  ) 


and 


m„a 


m,a 


m„a 


(nijSinh  m^x)  •  (cosh — 2 — )  -  (m2  cosh — j— )  *  (sinh — j-) 


a  =  Pmjm2at  • 


m.a 


m„a 


m.a 


m„a 


(m^  cosh — y~ )  *  (sinh — ^— )  -  (m2sinh — — )  •  (cosh — j-) 


where. 


m 


_  6E  . 

1  -  ~J~  (1  +  m3) 

i  yj 


m  =  -4^-  (1  -  m„) 
2  tG 


„  2G2t  v7 
m3  "  (l  '  SEEjtj  * 
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The  shear  stress  distribution  shows  no  appreciable  difference  from  that  obtained 
by  Volkersen’s  analysis,  except  that  at  the  two  ends  of  the  overlap  it  rapidly  falls  to 
zero.  This  of  course,  because  the  boundary  conditions  used  in  the  analysis  require 
that  the  shear  stress  be  zero  at  these  two  points.  Since  the  two  end  surfaces  of  the  ad¬ 
hesive  are  stress  free,  it  is  argued  that  for  equilibrium,  the  shear  stress  must  vanish 
at  the  two  ends.  However,  it  should  be  noted  that  the  lap  joint  has  a  discontinuity  at 
the  two  leading  corners  of  the  adhesive  and  the  stress  at  these  points  can  be  expected 
to  be  singular  for  an  elastic  analysis.  There  is  no  such  discontinuity  at  the  two  trailing 
comers  however,  emu  the  shear  stress  at  these  two  points  can  be  expected  to  be  zero. 
Since  Kelsey  and  Benson  do  not  allow  the  stresses  to  vary  over  the  adhesive  thickness, 
the  imposition  of  the  boundary  condition  that  the  shear  stress  be  zero  at  the  two  ends 
of  the  overlap  only  results  in  the  shear  stress  falling  rapidly  to  zero  at  these  two  points. 
The  shear  stress  distribution  cannot  however,  be  considered  an  improvement  over  that 
obtained  by  Volkersen. 

The  tearing  stress  distribution  obtained  by  Kelsey  and  Benson  shows  these 
stresses  to  be  concentrated  almost  entirely  at  the  two  ends  of  the  overlap.  Through 
most  of  the  length  of  the  overlap,  this  stress  is  almost  zero. 

2,6  Cornell’s  Modification  (10) 

The  system  considered  by  Cornell  is  shown  in  Fig.  4.  This  joint  is  fixed  at  one 
end  and  it  is  some  what  different  from  the  conventional  lap  joint.  Cornell’s  interest 
lay  in  brazed  joints  where  a  thin  tab  is  brazed  to  a  thick  base  bar.  The  base  bar  is 
subjected  to  tension,  bending  moment  and  shear  force.  Cornell  assumed  that  the  two 
adherends  act  like  simple  beams  and  the  adhesive  layer  consists  of  an  infinite  number 
of  shear  and  tension  springs.  Differential  equations  for  a  and  t  are  combined 
into  a  tenth  order  linear  differential  equation  in  the  transverse  deflection  of  each  ad- 
herend.  '  ’he  mathematical  analysis  is  fundamentally  simple  but  the  expressions  involved 
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are  complicated.  Results  indicate  a  high  concentration  of  shear  and  transverse  tensile 
stress  in  the  vicinity  of  the  joint  edge. 

Cornell  has  compared  his  theoretical  results  with  the  results  of  photoelastic  and 
brittle  lacquer  experiments  and  found  them  in  good  agreement.  It  should  however,  be 
noted  that  the  joints  investigated  by  Cornell  both  theoretically  and  experimentally,  are 
brazed  joints  where  the  adhesive  modulii  are  1/3  to  1/2  as  large  as  those  of  the  ad- 
herends,  whereas  in  an  adhesive  joint  the  ratio  might  be  1/20  or  smaller.  The  error 
in  the  predicted  stresses  due  to  the  assumption  that  the  stresses  do  not  vary  over  the 
chickness  of  the  adhesive,  may  be  more  serious  in  the  case  of  the  adhesive  joints  than 
it  would  be  for  brazed  joints. 

Cornell's  theory  applies  to  joints  between  adherends  differing  in  material  prop¬ 
erties  and  thickness.  Application  of  external  shear  R  :ce  and  bending  loads  are  also 
considered.  The  configuration  considered  by  Cornell  is  however,  of  a  special  type 
and  the  theory  cannot  be  applied  to  a  general  case  of  lap  joint. 

2.  7  Stepped  Joint  of  Erdogan  and  Ratwani  (11) 

In  this  analysis  a  stepped  joint  between  two  plates  is  considered  (Fig,  5), 

Plate  (1)  is  considered  to  be  isotropic  while  Plate  (2)  is  orthotropic.  A  uniform  tensile 
force  P  is  applied  to  the  plates  away  from  the  joint.  The  tensile  stresses,  Oj(x)  and 
a2(x),  in  the  two  plates  and  the  shear  stress,  x  ,  at  the  interface  are  determined  under 
the  following  assumptions: 

(1)  The  thicknesses  of  the  two  plates  are  small  compared  to  the  other  dimensions 

of  the  composite  structure  so  that  the  individual  layers  and  the  composite  plate  may  be 
considered  to  be  under  generalized  plane  stress  (i,e.  =  o2y  =  0). 

(2)  The  thickness  variation  of  the  stresses  in  the  plates  is  neglected. 

(3)  In  the  z  direction  (Fig.  5a),  it  is  assumed  that  either  c^z  =  e2z 


or 
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a  =0,6  being  the  average  stress  in  the  composite. 

z  z  r 

The  results  of  this  analysis  are  then  applied  to  a  special  case  of  one  single  step 

(Fig.  5b)  which  of  course  is  a  single  lap  joint.  The  shear  stress  in  the  adhesive  layer  is 

given  by, 

2  cosh  ax 

sinh  ax  +  (1  -  cosh  era  -h  -&-)  - ; - 

P  sinh  a  a 

where. 


a2  = 


1  - 


1  - 


^x  '  “2z 


Elrl 


E2xt2 


fi  = 


G(1  -  u{) 


For  both  plates  isotropic  and  identical. 


a 


2 


2 G 
t 


Eltl 


and 


T 


cosh  ax 

sinh  ax  -  (1  +  cosh  a  a)  - 

sinh  aa 


This  result  is  similar  to  that  obtained  by  Volkersen.  This  theory  assumes  that  the 
tearing  stresses  are  noneyistant  in  the  adhesive  layer. 
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2.  $  Finite  Element  Analysis 

Problems  in  structural  mechanics  can  be  numerically  solved  by  using  Else  finite 
element  method  (12).  Wooley  and  Carver  (13)  have  attempted  to  apply  tins  me  tired  to  Ebe 
case  of  a  single  lap  joim.  Tbe  lap  joint  is  loaded  tn  tension  and  stesr  and  tearing 
stresses  are  determined,  Tire  adhesive  layer  was  divided  into  two  equal  layers  through 
its  thickness,  giving  an  indication  of  the  variation  of  stresses  in  me  direction  normal 
to  the  bondline.  The  solution  is  obtained  by  a  plane  stress  analysis.  Tire  shear  and 
tearing  stress  distribution  obtained  by  this  method  is  similar  to  that  obtained  by  Golanri 
and  Reissner,  Wooley  and  Carver  compare  the  stress  concentration  factors  wim  those 
obtained  by  Golsnd  and  Reissner  theory.  They  are  shown  to  compare  favorably. 

The  joim  considered  by  Wooley  and  Carver  is  loaded  in  tension  only.  The  finite 
element  method  is  however,  a  promising  one  and  can  be  extended  to  a  more  general  case 

Another  attempt  to  use  the  finite  element  analysis  in  the  case  of  adhesive  joints 
is  tha‘  by  Harrison  and  Harrison  (14).  It  is  assumed  that  the  displacements  at  the  two 
interfaces  are  known  and  are  uniform.  This  is  a  very  special  case  and  the  closest 
practical  case  would  be  one  where  the  adherends  are  rigid  It  is  funner  assumed  that 
in  a  region  of  the  adhesive  layer,  away  from  the  two  ends,  the  stress  field  is  uniform 
This  uniforn  stress  field  is  determined  from  the  known  displacement  field  at  the  two 
interfaces  f-nite  element  analysis  is  then  applied  to  determine  stresses  near  the  ends 
of  the  adhesive  layer. 

In  practice,  the  relative  displacements  of  the  interfaces  are  not  known,  only 
the  external  loads  rre  given.  Furthermore,  the  displacements  at  the  interfaces  are 
far  from  being  uniform,  especially  at  the  ends  of  the  joint  This  method  needs  to  be 
further  developed  in  order  to  be  applicable  to  practical  problems  in  adhesive  joints. 
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2. 9  Experimental  Stress  Analysis 

Several  works  have  been  published  where  the  stress  distribution  was  determined 
experimentally.  These  works  will  be  reviewed  here  briefly. 

Unpublished  extensometer  tests  made  by  Copper  have  been  referred  to  by  Mylonas 
and  deBruyne  (5).  In  these  tests,  the  shear  stress  developing  at  various  points  of  the  3 
inch  long  overlap  was  estimateo  from  the  relative  displacement  of  the  adherends  meas¬ 
ured  by  an  optical  lever.  The  load  applied  was  tensile  and  was  kept  very  low  so  that 
Hooke's  law  could  be  assumed  to  hold  even  in  the  region  of  high  stress.  Tne  stress 
distribution  so  obtained  was  compared  with  the  results  of  Volkersen's  theory  and  also 
with  the  results  of  Goland  and  Reissner  theory.  Experimental  results  were  found  to  be 
in  better  agreement  with  the  theory  of  Goland  and  Reissner  than  with  the  theory  of 
Volkersen. 

Mylonas  (7)  applied  the  method  of  photoelastic  analysis  to  investigate  the  stresses 
in  a  lap  joint.  The  adherends  were  made  of  \  in.  square  Bakelite  bars  reinforced 
longitudinally  with  stainless  steel  studdings.  The  adherends  were  made  practically 
rigid.  The  adhesive  used  was  an  Araldite  resin  which  is  also  a  photoelastic  material. 

The  thickness  of  the  adhesive  layer  was  appreciable. 

The  main  purpose  of  Mylonas'  work  was  to  study  the  effect  on  the  peak  stress 
of  the  shape  of  the  free  boundary  of  the  adhesive  at  the  edge  of  the  overlap.  Models 
with  concave  edges  with  varying  radii  of  curvature  and  with  straight  edges  of  various 
inclinations  were  studied. 

The  practically  infinite  stiffness  of  the  adherends  in  these  experiments  leads 
theoretically  to  a  uniform  shear  throughout  the  adhesive  layer.  For  this  reason  Mylonas 
results  are  of  limited  interest  from  the  stress  analysis  point  of  view.  The  work  is 
however,  of  great  importance  as  a  pioneering  work  describing  techniques  for  making 
stress  free  adhesive  joints  to  be  used  as  models  for  photoelastic  analysis. 
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McLaren  and  Maclnnes  (15)  studied  the  effect  of  bending  of  the  adherends  in  their 
photoelastic  tests.  Two  series  of  tests  were  made..  In  the  first,  a  lap  joint  consisting 
of  5  in.  thick  adherends  and  an  adhesive  layer  varying  in  thickness  from  zero  to  \  in. 
was  simulated  by  a  model  cast  whole  in  Araldite;  in  the  second  series  of  tests  the  ad¬ 
herends  had  a  Young's  modulus  twenty  times  greater  than  that  of  the  adhesive  layer. 

In  the  homogeneous  models  the  'adhesive  layer'  was  always  rectangular  and 
ran  the  whole  length  of  the  overlap.  In  every  case  in  which  k  >  0  the  fringe  pattern 
had  the  same  general  character:  the  order  was  a  minimum  at  the  center  of  the  joint 
and  increased  towards  a  maximum  at  the  leading  comer.  Change  of  the  value  k  pro¬ 
duced  an  approximately  proportional  change  in  the  value  of  the  peak  fringe  without 
appearing  to  affect  the  order  at  the  center.  For  negative  (k  =  -  \),  the  fringe  order 
at  the  free  end  of  the  glue  line  reduced  to  a  minimum,  the  fringe  order  was  a  maximum 
at  the  center. 

In  composite  models  which  were  geometrically  similar  to  the  homogeneous  ones, 
similar  fringe  patterns  were  obtained.  For  k  >  0  the  fringe  order  rose  from  a 
minimum  at  the  center  of  the  adhesive  layer  to  a  maximum  at  the  leading  corner.  The 
high  tensile  stress  in  the  free  surface  at  the  leading  corner  decreased  on  traversing 
the  end  but  gave  way  to  slight  compression  at  the  other  extreme.  Increase  in  the  value 
of  k  again  produced  no  change  at  the  center  but  led  to  an  increased  fringe  order  at 
the  leading  comer.  When  k  was  made  sufficiently  negative,  it  was  found  that  the 
highest-order  fringe  appeared  at  the  center  of  the  adhesive  and  then  there  was  a  fall 
in  order  towards  the  free  surfaces. 

The  main  significance  of  this  work  by  McLaren  and  Maclnnes  is  that  it  clarifies 
the  contribution  made  to  the  distribution  of  stress  by  the  bending  of  adherends.  The 
results  tend  to  support  the  Goland  and  Reissner  theory, 

Tuzi  and  Shimada  (8)  use  both  the  adhesive  and  the  adherends  made  of  photo- 
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elastic  materials.  The  adherends  are  made  of  an  epoxy  resin  while  the  adhesive  is 
made  of  epoxy  rubber.  The  ratio  of  the  Young’s  modulii  for  the  two  materials  is  6. 

The  fringe  pattern  obtained  for  rectangular  adhesive  layer  is  similar  to  that  obtained 
by  Melaren  and  Machines.  Tuzi  and  Shimada  study  adhesive  layers  with  concave  and 
canvex  ends  and  joints  with  a  fillet  at  the  leading  comers.  The  results  indicate  that 
convex  ends  tend  to  produce  excessive  stress  concentration  at  the  leading  comer  where¬ 
as  the  use  of  a  fillet  tends  to  reduce  the  stress  concentration  considerably. 

The  photoelastic  investigations  by  Mylonas,  McLaren  and  Machines,  and  Tuzi 
and  Shimada,  all  show  a  significant  variation  of  stress  through  the  thickness  of  the 
adhesive  layer,  especially  in  areas  where  the  maximum  stress  occurs. 

2.10  Conclusion 

Amongst  the  theoretical  works,  Goland  and  Reissner's  theory  is  the  most 
rigorous  study  of  the  stress  distribution  in  lap  joints  and  is  in  good  agreement  with 
experimental  results  except  at  and  near  the  two  ends  of  the  joint  where  the  maximum 
stress  occurs.  This  error  results  mainly  from  the  assumption  that  the  stress  does 
not  vary  through  the  thickness  of  the  adhesive  layer.  In  fact,  photoelastic  tests  show 
this  variation  to  be  quite  significant  especially  in  those  areas  where  the  maximum 
stress  occurs.  Furthermore,  the  theoretical  analyses  consider  the  joint  to  be  loaded 
in  tension  alone  and  the  effects  of  shear  force  and  bending  moment  are  not  studied. 
Photoelastic  stress  analysis  by  McLaren  and  Machines  shows  that  application  of  shear 
force  and  bending  moment  can  be  advantageous  from  the  point  of  view  of  obtaining  a 
more  uniform  stress  distribution  along  the  overlap. 

The  objective  of  the  present  work  is  to  analyze  the  lap  joint  subjected  to  a 
generalized  load  configuration  consisting  of  tensile  and  shear  forces  and  bending 
moments  (Fig.  3).  The  stresses  are  allowed  to  vary  through  the  thickness  of  the  ad¬ 
hesive  layer.  The  analytical  results  are  to  be  checked  by  means  of  photoelastic  tests. 
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3.  ANALYSIS  OF  LAP  JOINT 


3, 1  Problem 

Given  are  two  rectangular  sheets  of  thickness  t^  and  t£  and  of  unit  width,  lap- 
jointed  over  a  length  'a'  by  means  of  an  adhesive  layer  of  thickness  t  (Fig,  So).  It 
is  assumed  that  the  loads  P,  V,  ,  and  M2  are  given  at  the  joint  edges.  The  main 
problem  is  the  determination  of  the  shear  stress  r  ,  and  the  normal  stresses  ax 
and  a  in  the  adhesive  layer. 


3.2  Assumptions 

The  following  assumptions  are  made: 

(1)  The  adherend  and  the  adhesive  materials  are  isotropic  and  linearly  elastic. 

(2)  The  adherends  behave  according  to  simple  beam  theory,  i.e.  the  lateral 
stress  in  the  adherends  can  be  neglected,  and  plane  sections  remain  plane, 

(3)  The  adhesive  thickness  t  is  small  compared  to  the  overlap  length  a. 

(4)  The  displacements  in  the  joint  are  small,  i.e,  the  strains  are  considerably 
less  than  1, 

(5)  The  joint  width  is  large  compared  to  adherend  and  adhesive  thicknesses, 
so  that  this  may  be  considered  as  a  problem  in  plane  strain. 


3.  3  Derivation  of  Differential  Equations 

The  coordinate  axes  x  and  y  are  shown  in  Fig.  3.  The  x~axis  is  parallel  to 
the  plane  of  the  adhesive  layer  <_nd  passes  through  its  center.  The  y-axis  is  normal 
to  the  plane  of  the  adhesive  layer  at  the  left  edge  of  the  joint. 

Let  u(x,  y)  and  v(x,  y)  represent  the  displacements  of  the  point  (x,  y)  in  the 
adhesive,  in  the  x  and  y  direction  respectively.  Further,  let 


u(x,y)  = 


u  (x) 
0  '  ' 


+  yUj  (x)  + 
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and 


v(x,y)  =  vq(x) 


yvj(x) 


+  y  v0  (x)  + 
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2 
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(3,2) 


where  m  and  n  are  integers. 

A  perfect  solution  would  require  m  and  n  to  approach  infinity.  However,  since 
the  adhesive  thickness  t  (dimension  in  the  y-direction)  is  small  compared  to  the  over¬ 
lap  length  a  (dimension  in  x -direction),  m  and  n  could  be  terminated  at  values  2,  3, 
or  4  and  the  results  could  be  expected  to  be  reasonably  accurate. 

The  strains  in  the  adhesive  layer  can  new  be  written  as. 


e 


x 


9  u  (x,  y) 
Sx 


m 

2  71  u.  (x) 

j  =  0  J 


(3.3) 


=  9  v(x.y) 

y  ay 


2  jy^'^v.(x) 

j  =  i  J 


(3,4) 


9  u(x,  y)  ,  a  v(x,y) 
3y  3  x 


11  ('  X 

2  jyu  V(x)  +  2 

j  =  1  J  j  =  o 


yj  v .  (x) 


(3.5) 


For  the  plane  strain  problem,  the  strain  energy  I  of  the  adhesive  is  given  by. 
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Integration  over  the  thickness  of  the  adhesive  yields. 
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The  condition  of  compatibility  at  the  two  interfaces  can  be  used  to  determine  the 
displacement  field  and  the  strain  energy  of  the  two  adherends. 

Let  Uj(x,y)  and  (x,y)  represent  the  displacements  in  the  x  and  y  direc¬ 
tions  respectively  in  adherend  (I). 

Since  plane  sections  remain  plane, 

9  V. 

Uj(x,y)  =  UL  (x,  -f)  -  (y  -  y )  (3.8) 

Now,  conditions  for  compatibility  at  the  interface  are. 


U^x,-!)  =  u  (x,  ^-  )  (3.9) 


and 

Vj  (x,  y )  =  v  (x,  y  )  (3.10) 


Eqs.  (3.  8),  (3. 9)  and  (3. 10)  yield 

_t_ 

Uj  (x, y)  =  u(x,-L)  -  (  y  *  y  )  --V-  (^-2  }  (3.11) 
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Similarly  for  adherend  (II)  the  displacements  U2  (x,y)  and  V2  (x,y)  in  the  x  and  y 
directions  respectively  can  be  written  as, 
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The  strain  field  for  adherend  (I)  is. 
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j  =  0  3 


(3.14) 


y.  =0 
lxy 


(3.15) 


And  the  strain  field  for  adherend  (II)  is. 


tn  »  t 

e2x  =  ,sn  <'T))u'i  '  fr  +  T> 

J  =  0  J 


2  (  -  4-  )3V". 

j  =  0  3 


(3.16) 


'y2xy  =  0 


(3.17) 


26 


The  strain  energy  1^  for  adherend  (I)  is  given  by, 


*  (4-+v  o 

■,-j  v  i 


e^  dy  dx 


4  <Tl  V 


T 

j  j  i 


O  -7T 


m 


2  (4-)Ju'.  -  <y-4-) 


n 


L  j  =  0 


2  j 


s  (-T>jv"j 


j  =  0 


dy  dx 


(  2  i  V  ql 

J  ~ 


_t_ 

2 


m 

2 


m 


2  (-§-)i+ju\u\ 

i  =  0  j  =  0  1  J 


+  (y  -  4->2 


n 

2 


n 


i  =  0  j  =  0 


s  (^-)i  +  jv"iv"j 


(~  m  n  .  . 

2(y4)  2  2  +  J  U’v" 

2  ( i=0  j=0  2  J  JJ 


dy  dx 


Integrating  in  the  y -direction. 


m 

2 


m 


i  =  0  j  =  0  1  3 


+ 


n 

2 


i=0  j  =  0 


2  (-y)i  +  j  vVj 


m  n  .  . 

2  2  (4)1  +  J  u’  v". 

i  =  0  j  =  0  z  1  J 


dx 


(3, 


13) 
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Similarly,  the  strain  energy  I2  for  adherend  (II)  is  given  by, 

a  Q2  f  m  m  i+i 

t  _  C  v  v  «-  „•  n* 


f  - 

Jo  2 


o 

m 

m 

l 

) 

2 

2 

1! 

O 

j  =  o 

n 

2 

n 

2 

(“T 

i  = 

0  j  =  0 

m 

n 

2 

2  ( 

-J-)1 ' 

=  0 

j  =  o 

(3.19) 


Work  done  by  the  boundary  forces  P  and  V  and  moments  and  can  now  be  calcu¬ 
lated.  Let  Wpl ,  Wy,  ,  and  represent  the  work  done  by  P,  V,  and  respectively 
on  adherend  (I).  And  let  Wp2 ,  Wy2  ,  and  WM2  represent  the  work  done  by  P,  V, 
and  M2  respectively  on  adherend  (II). 

Assuming  P  to  be  uniformly  distributed  over  the  cross-sectional  area  of  each 
adherend, 


<T  tV 


—  U,  (a,y)dy 
z\  1 


(T  j  tl) 


m  (  n  .  ) 

S  (4-)Ju.(a)  -  (y --f)  ]  2  ( L  dy 

j  =  o  2  ( j  =  o  J  r 


m  t  i  ci  n  t  i 

P  2  (4-)3u(a)  -  -±-  2  (4)Jv’.(a) 

j  =  0  2  ]  2  j  =  0  2  J 


(3,20) 
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Similarly, 


W  =  -P 
P2 


m 


j  =  0 


*2  n 


S  (--l)Ju(o)  +  2  (-f>Jv\-<°> 

2  J  2  j  =  0  Z  J 


(3.21) 


W\<1  =  V'Vl<a)  =  V  2  (4-)jv?(a) 


j  =  0 


2'  j' 


(3.22) 


n 


WV2  =  -V*  V2(o)  = 


-V  2  (-J-)V(o) 


j  =  0 


2  7  j 


(3.23) 


wM1  =  M,  dvJ*’J/21 

Ml  I  3x 


n 

=  M  2  (-L)V(a) 
1  j  =  0  Z  J 


(3.24) 


WM2  =  -M*>  =  (-4^*^) 

M  “  j=0  “  J 


(3.25) 


The  total  work  done  by  external  forces  is  given  by, 


W  =  W  +  W  +  W  +  W  +  W  +  W 

PI  P2  VI  V2  Ml  M2 


(3.26) 


Let  Iq,  Ip  and  ^  be  expressed  as. 


In  =  j  Fn(u. ,  u\  ,  v. ,  v*  ,  v". )  dx 
0  J  O'  i’  i»  j'  j’  j ' 


(3.27) 


d 

I,  =  f  F.  (u. ,  u’. ,  v. ,  v’. ,  v". )  dx 
1  Jn  r  i*  l  J  j’  V 


(3,28) 


L  =  i  F0  (u. ,  u'. ,  v. ,  v\ ,  v". )  dx 

2  J  2 '  T  i’  j’  j’  j ’ 


(3.29) 


where  i  =  0,  1,  2, . .  m  and  j  =  0,  1,  2,  n 


i 


The  functions  F  q  ,  F^ ,  F2  are  given  by  the  expressions  on  the  right  hand  side  of  the 


equations  (3. 6),  (3. 18),  and  (3. 19)  respectively.  Let 


F(u. ,  u'.,  v.,  v'.,  v".)  =  Fn  +  F.  +  F0 
v  1  1  j  j  j  0  1  2 


(3.30) 


so  that  the  total  potential  energy  1  can  be  written  as. 


f  F(u.,  u\,  v.,  v'.,  v".)  dx 
Jn  1  1  J  J  J 


(3.31) 


Now,  equilibrium  requires  that  I  be  minimized.  Application  of  the  principle  of  calculus 
of  variations  yields  the  following  differential  equations;  also  known  as  Euler  equations. 


—  (— I— )  =  0 

dx  v  au'.  ’ 

1 


i  =  0,  1,  2,  ,  m 


(3.32) 


9F  d  ,  9F  ,  ,  d  /  9F  N  _  n  .  _  n  ,  ,0 

9v.  "  'dx  ^  9v*.  *  +  ~7T  ^  9v".  ^  0  J  0>  l’  - - -  •  n  ^3*33) 

J  J  dx  0  j 


In  order  to  eliminate  the  rigid  body  displacements  of  the  joint,  it  will  be  convenient  to 
let  the  displacement  and  slope  at  the  origin  lie  -  ero.  This  would  yield  the  forced 


boundary  conditions  (3. 34). 


uQ (0)  =  0 


vQ(0)  =  0 


v’p(0)  =  0  (3.34) 


The  natural  boundary  conditions  are  given  by  Eqs.  (3.35)  to  (3.40). 


-P(-T>i  =  0 


i  1,  2,  »...o.., 


m 


(3,35) 
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(— dF  _)  _  p(-L)1  =  o 

'  9u\  ;x  =  a  '  2}  v 
1 


i  =  0,  I,  2, 


(3.36) 


'  av\  x  =  o 
} 


dx  K  dw".  1 


x  =0 


-  V(-4)j  =  0 


j  =  1,  2, 


(3.37) 


'  av’j  ;x  =  a 


JL  (_9£_) 

dx  V  av’\  1 


-  Vi-k-?  = 


j  =  0,  1,  2, 


(3.38) 


h$->=c  =  0  -^-¥  (-T-  +  M2>  -  0 


j  =  1,  2,  . . . . ,  n  (3. 39) 


(— 9_ ...)  +  (— )^( — — 

1  9v".  ;x  =  a  ^  2  M  2 
J 


-  Mj )  =  0 


j  =  0,  1,  2 . n  (3.40) 


Differentiation  of  F  yields  Eqs.  (3. 41)  to  (3.  50). 


(3.41) 


■fr  -  G  iJ<V+')'1  6ii  'j  +  '.V<t^+1  vj 

l  j=l  j=0 


i  =  1,  2, 


(3,42) 


_3F  _  Q_ 

3U1-  “  2 
l 


S  2(-L)i  +  j  +  I  «  u- 

j  =  0  ~  1J  j 


+  A  2  j( 4- )1+j  /3..V. 

L  j  =  o 
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Qr 


E2S  .  ^  a 

[  j  =0  1  "  J  B  j  =0  -  J 


02 


EB 


z  2z,(~r  3u’ 
j  =  0  -  - 


*  s“ 


Q 

y  l  — 7  J 


1  = 


i  =  0,  1,  2, 


. ,  ins 


(3.43) 


—  (  dF.  ) 


=  4  I  .r  <-?-> 

j  =  0 


HI 


Z.1+} 


a . .  u. 

*J  J 


Xt 

9 


r  i  —  i  -  1 


S  j(J_)—  J  *  0..v*. 

j  =  i  2  ,J  Jj 


.  Qi'i 


^2r2 


m 


S  2(4-/  '*  j  u”  -  t, 
j  =  0  Z  J 


n 

S 


i =  o 


<4 :  j  '■'"j 


m 


2  2  (-4) 

j  =  0 


t  a  +  j 


u  . 
J 


t,  £  (-4)i:jv- 

”  j  =  0  -  J 


i  =  0,  1,  2, 


,  m 


_9F 

dv 


-  =  0 


aF  _ 


a- 


Q_ 

2 


(3.44) 

(3.45) 


n  ... 

1 

~ 

Z  2ij(i)1  +  )  "  6  V. 

1=1 

+  \ 

i  =  0  J  J 

j  =  1,  2, 


(3.46) 
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n  -  n 

j  =  l  2  03  3  j=0 


(-r)t22('1)j 


Ql  2 

-  <4->y 


<T>iv’"j  =  0 


(3. 51) 


m  t  i  -  1  m 

G  S  ij(-W  6  u  -  £ 

1=1  2  1J  J  J-0 


^r“.j  +  QA  +  Q2t2(-Di  + j 


n 


n 

2 


+  S  (Gi -Xj)(-f)]/3  V 

j  =  0  2  13  3  j  =  o 


<-?)t22(-l)i+3 


.  (— L)t 2  ■  (— )3  v’"  =  o 
'  2  '  1  j  '  2  '  v  j 


i  =  1,  2, . .  m 


(3.  52) 


Differential  Eqs.  (3.33)  can  be  rewritten  as  Eqs.  (3. 53)  and  (3.  54). 


m 


G  2  i(4-)ii3i  u'.  +  2 
i  =  1  2  10  1  i  =  0 


m 


<^->.22  (-D1  -  4->'i2  i 


■G  S  (-§•>*  +  1  “io  v"i  +  = 
i  =  0  2  10  1  i  =  0 


Ql'l3  Q2t23  i 

~y—  +  (-D 


,  t  si  iv  n 

<-)v.  -o 


(3.53) 


m 


2  (Xj  -  GiX-^-)1  /3  u'.  +  S 


m 


i  =  0 


i  =  0 


.  ^2..  2.  ,  ,i+  i  A  ,,  2 

(  — )t2  (-1)  -  (—  )t, 


<T>‘u'"i  +  <f>  .  =  2‘j  <T>‘  “  6i  j'1!  -  0 


Q 


t  \i  +  1 


n 

+  2 
i  =  0 


Q^l  ^2^2  i  +  i 

1 1  +  -y~-  (-i)1  +  J 


(T)iyiiV  =  0 


j  =  1,  2, 


(3,  54) 
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The  forced  boundary  conditions  are, 


u  (0)  =  0 
o  v 


v  (0)  =  0 
o  ' 


v*  (0)  =  0 
o 


The  natural  boundary  conditions  can  be  rewritten  as  Eqs.  (3. 56)  to  (3. 61). 


m 

2 

j  =  0 


n 


<i>V  (0)  +  s  xj<-§.)J0  v 

.1  J  =  1  J  J 


n 

+  2 
j=0 


(^-)t22(-l)i  +  )  -  (-J-)!; 


(■y)jv".( 0)  -  (-l)P  =  0 


i  1,  2,  ......a,  m 


m 

2 

j  =  0 


(  2  )ta  +  Q1t1  +  Q2t2  (-1) 


i  +  j 


n 


(-f)J  u'.  (a)  +  2 

J  J  =  1 


Xj  (-i-)j^  (a)  +  2 

J  j  =  0 


^2  2  i  +  i  ^1  2 

{~T)l2  (~l)  J  -  ("2L)tl 


(i>j  vya) 


i  “  Oj  Ij  2  f  »»«|«Q03f  1T1 


m  i  m 

2  (3.  ■  u  (0)  -  2 

i  =  1  Z  13  1  i  =  0 


.  Q2X  2  .  i  +  j  Q1  .2 
(  2  ^2  ^  (  2  ^1 


(±.fu"(0)+  2  C(4)i  +  1  a  v1.  (0)  -  2 

L  i  =  0  Z  1J  1  i  =  0 


Qiti3  ,  V23 
— _ —  +  — * —  ^ 


(t^’V0*  "  ('1)v  =  0 


j  ~  1,  2,  o.o.ooo,  n 


(3.55) 

(0) 

(3. 56) 

-P  =  0 

(3.  57) 


(3, 58) 
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m  .  m 

S  Gi(-f)1j3.iu  (a)  -  2 

i  =  l  2  1J  1  i  =  0 


Qo  9  i  +  i  ^1  2 


+  2  G(-^-)i  +  1  a.  .v'.(a)  -  2 

i  =  0  J 


i=  0 


(^-)V':(a)-V  =  0 


j  =  0,  1,  2, 


,  n 


(3. 59) 


m 

2 

i  =  0 


,Q2  2  .  ..i+j  ,Q1  w2 
<— >t2  <_1)  <~)tl 


<T^u’i<0> 


n 

+  2 
i  =  0 


Qit,3  Q2t23  i  +  i  t  i  ;  Pt2 

~Y~  +  +  J  (-|-)1v"i(0)  -  (-1)J(-/-  +  M2)  =  0 


j  —  1,  2,  . .  n 


(3.60) 


m 

2 

i  =  0 


,^2  2  .  ,,i  +  j  ,  ^1  .,2  ,  t  ,i ...  ,  i 

(~)'2  (-D  -  >ll  <T)ui(a) 


n 

+  2 
i  =0 


<T>‘v"i<a>  +  <"£• 


-  Mj)  =  0 


j  =  0,  1,  2, 


,  n 


(3.61) 


3.4  Parameters  of  Similarity 

Inspection  of  the  differential  equations  and  the  boundary  conditions  shows  that 
the  displacement  field  in  the  adhesive  layer  is  a  function  of  G,  X,  Q^,  Q2>  t,  tj,  t2,  a, 
P,  V,  M^and  M2<  The  equilibrium  of  the  external  forces  however,  gives  M?  as  a 
function  of  P,  V  and 


*1  ^2 

M2  =  Mj  +  Va  -  P(-j-  +  -g-  +  t) 


(3.62) 
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Further,  it  can  safely  be  assumed  that  the  Poisson's  ratio  4  does  not  vary  appreciably 

for  different  adhesive  materials.  G  and  X  could  then,  be  considered  as  functions  of 

Q.  Also,  the  various  dimensions  and  the  loads  can  be  expressed  in  dimensionless  form. 

t  *1  *2  P  V  ^1 

The  parameters  of  similarity  then,  are,  Q,  Q,  »  Q0  ,  —  ,  - ,  - ,  — ,  -p-,  and  -p—  » 

L  Z  3.  3.  cL  3.  x  a  3, 

The  stresses  a  and  a  are  dependent  on  the  modulus  of  elasticity  E,  while 

x  y 

the  shear  stress  t  is  determined  by  the  shear  modulus  G.  Since  4  does  not  vary 

greatly,  G  can  be  considered  as  a  function  of  E.  Also,  if  we  assume  that  4,4,,  and 

Q1  Q1  E1  Q2  Q2  E2 

42  are  approximately  equal  then,  -g-  a  -q-&  -g-  and  -g-  a  -g--a-g- .  The  stress 

^1  t  t. 

field  in  the  adhesive  layer  is  then  determined  by  the  factors,  — p-,  — =r- ,  —  ,  — —  , 
t  „  M  b  b  a  a 

a  *  a  *  P  md  Pa  * 

A  special  case  of  identical  adherends  with  antisymmetric  loading  configuration 
can  be  considered  here.  In  this  case, 


Mj  =  -M2  =  \ 


P(t1  +  t)  -  Va 


The  stress  field  in  this  case,  is  a  function  of 


and 

!l  j 

'E  ’  t, 


h  l2  ’  E1  E2 


1  P  V 
a  ’  a  ’  P 


(3.63) 


3.  5  Case  of  Identical  Adherends 

The  differential  equations  and  the  boundary  conditions  can  be  greatly  simplified 
if  the  two  adherends  are  identical  and  m  and  n  are  taken  to  be  equal  to  1.  An  analytical 
solution  of  the  differential  equations  is  possible  in  this  case.  Since,  the  adherends  are 
identical, 


t2  =  tj  ;  and  Q2  =  Q[ 


(3.64) 


The  four  differential  equations  in  this  case  are: 


(Qt  +  2  Q[t1  )uq"  +  Xtvj' 


tv 


1 


1 1  1 


0 


(3.65) 
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2GUJ  -  <^-  +  2Q1t1)^-  u,”  +  2Gvo'  +  QftV  =  0 


(3.66) 


GnV  -  ^1*1  T  “l  ”  -  GtV  +  ^1  ~5  voiV  '  0 


(3.67) 


2  t,3t 

aV  -  QiV “o'"  +  ^  -  ttV  +  vi*v  = 0 


(3.68) 


The  boundary  conditions  in  this  case  are  given  by  Eqs.  (3. 69)  to  (3. 78). 


u  (0)  =  0 
o  ' 


vo(0)  =  0 


vo’(0)  =  0  (3.69) 


(%-  +  2Q1t1)-fu1,(0)  '  Qi*iV<0)  +  P  =  0 


(3. 70) 


t.2t 


(Qt  +  2Qit1)uo’  (a)  +  Atvj(a)  -  QL  — 2— V  (a)  -  P  = 


(3.71) 


(-S1  +  2Q1t1)-i-u1'(a)  -  Q1t12v0"  (a)  -  P  =  0 


(3.  72) 


Gt„ 


tL  t 


Qltl2uo"(0)  +  -^2.Vl'  (0)  -  Q1-ir_v1"(0)  +  V  =  0 


(3.73) 


t2t 


2t, 


GtUj  (a)  +  Q1_l_u1"(a)  +  GtVQ' (a)  -  Q1_^_vo"'(a)  -  V  =  0  (3.74) 


2  t  2t 

Q^V'ta)  +  -T“V(a)  "  3l-V-vl”'(a>  '  v  =  0 


(3.75) 


t,3t 


Pt, 


~C*ltluo  (0)  +  Ql-V-Vl"(0)  +  <-T"+  M2>  = 


2 


0 


(3.76) 


38 


t,2t 


t  2t  Pt 

•Ql-T-'V<a>  +  liT'.'W  +  <-/  -  Mj)  =  0 


(3„  77) 


tA 


Pt, 


"Q^Uq’ (a)  +  Qi-! — Vj"  (a)  +  (-j-  -  Mj)  =  0 


(3. 78) 


The  solution  for  the  differential  equations  can  be  found  as  follows: 
From  (3. 65)  we  have, 


uo"  =  Aj  •  vf*  +  A2  •  vj* 


(3, 79) 


where, 


A  _  A*1  1 

A1  =  2(0t  +  2  ) 


and 


A 


2 


-\t 

(Qt  +  2Q1t1) 


Integrating  (3.  79), 


Uo’  =  A1  *  Vl"  +  A2  *  V1  +  C1 


(3.  80) 


(3.  81) 


(3.  82) 


where  is  a  constant. 

Substituting  (3.  79)  and  (3.  82)  in  (3. 68)  and  rearranging  the  coefficients,  we 
get, 


+  Vl" 


+  A4Vj 


c 


1 


(3.  83) 
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where. 


Gf 


A3  "  (~6 - 4XV  7  A5 


A4  =  -2(XA2  +  Q)  /  A5 


A5  AlQlfl 


2 


Qili 1 


(3. 83)  is  an  ordinary  linear  differential  equation  with  constant  coefficients  and  yields, 


k^x  -kjX  k2x  -k2x  2  X  c^ 

V1  =  c2e  +  c3e  +  c4e  +  c5e  +  A4A5 


(3.  84) 


4  2 

where  k^  ,  -k^ ,  k2  ,  -k2  are  roots  of  the  polynomial  k  +  A^k  +  A4  =  0  ;  and 

c2  ’  C3  ’  c4  ’  c5  are  constants  integration. 

Substituting  (3.  84)  in  (3. 82)  we  get. 


2  klx 


2  ‘klx 


2  k2X 


2  "k2X 


uo’  =  Al<c2kl  e  +  c3kl  e  +  c4k2  e  +  c5k2  e  > 


k.x  -k.x  k„x  -k„x  2Xc. 

+  A2  ( c2e  +  c3e  +  c4e  +  c&e  +  j-g-  )  +  ^ 


or 


A2C2  klX  A2C3  “kix 

uo  =  <V2ki  +  >e  '  <V3ki  +  -rf  > 


A9c„  k.x 

+  <Alc4k2  + 


A9c,-  -k9x  2A9X 

( Alc5k2  +  ~kj~'^e  +  (A^T  +  1)cix  +  c6 


(3.  85) 
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where  is  a  constant  of  integration. 

Now,  and  uq  are  known  functions  of  x  from  (3.  84)  and  (3.  85)  if  the  constants 
ci  *  c2  *  c3  *  c4  >  c5  *  anc*  c^are  known.  These  constants  can  be  determined  from  the 
boundary  conditions  (3. 69),  (3.71),  (3.73),  (3.75),  (3.76),  and  (3. 78),  These  boundary 
conditions  yield  the  following  set  of  simultaneous  equations. 

From  (3. 69),  i.  e.  for  uq  (0)  =  0, 

^Alkl  +  )C2  "  (Alkl  +  kj"  )c3  +  (Alk2  +  kj  )c4  "  (Alk2  +  kj  *C5  +  c6  =  0 

(3.  86) 

From  (3, 71), 

(Qt  +  2Q1t1)c1  =  P  (3.87) 

From  (3. 73), 

(A6ki  +  A7kj3)  c2  '  (A6ki  +  A7ki3)c3  +  (Agk2  +  A7k23>c4  '  (A6k2  +  A7k23^c5  =  'V 

(3.  88) 


where, 


A,  = 


¥- +  w 2 


and 


2  ^ 

Ay  =  Q1^  (Al  "  -5—  ) 
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From  (3. 75), 

3  ^1a  3 

(Agkj  +  A^  )e  1  c2  -  (A6k1  +  A^  )e  Cg 

3  ^?a  3  ^2a 

+  (A6k2  +  A7k2d)e  2  c4  -  (A6k2  +  A^e  2  c5  -  V  (3.89) 

From  (3. 76), 

+  (A^2  +  Ag)c2  +  (Aykj2  +  Ag)Cg  +  (A^2  +  Ag)c4 
2  Ptl 

+  (Ayk22  +  Ag)c5  =  (-y-+  M2)  (3.90) 

where, 

A8  =  _<Vl  A2 

and 

2AA„  9 

A9  =  -<1+  -ApC-^lV 

4  5 

From  (3,  78), 

2  ^la  2  _k.a 

A9c1+^L?k1  +  Ag)e  c2  +  (A^j  +  Ag)e  c3 

2  k_a  2  -k?a  Pt. 

+  (A?k2z  +  Ag)e  2  c4  +  (Ayk22  +  Ag)e  c5  =  Mi  "  -J~  (3.91; 

The  above  six  simultaneous  linear  equations  can  be  solved  to  obtain  the  six  constants  which 
when  substituted  in  (3,  84)  and  (3, 85)  would  give  v^(x)  and  uq(x)0 
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&ow,  to  determine  Uj(x)  and  vq(x),  differential  equations  (3. 66)  and  (3, 67) 
sfescsM  fee  solved. 

Differentiating  (3. 66)  once,  multiplying  by  t/2  and  then  adding  to  (3. 67)  gives, 

•Vfhr  •  i>voiv  -  |  ^  +  2QiV  4- +  V K-  V  = 0 


or 


S'  =  Vo 


r? 


(3.92) 


where. 


Oi 

9  -t. 


Bi  = 


*Vl  (  3  +  2  * 


<TT  ♦  2Q,V -T-  +  Q, 


t  1 


Integrating  (3. 92),  we  have, 

"l  =  Blvo'  +  ~T~  *2  +  V  +  CS 


(3.93) 


where  c_  ,  Cg  ,  are  constants  of  integration. 

Substituting  (3. 93)  in  (3. 66)  and  rearranging  coefficients,  we  get. 


v  +  B0v  '  =  B0c_x”  +  2B~c0x  +  2B0cn  +  B.c., 
o  2  o  37  38  39  47 


(3.94) 


where. 


B2  B 


2G(Bj  + 1) 
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B4  ’  - 


(ir-  +  Ofa* 

Br 


B5  ^T-+(^ltPtBl  "  ^Vl 


The  solution  of  (3. 94)  is. 


B3  3  B3  2  \  .  B4 

Vo  C7  3B0  X  +  C3  B0  +  1  (  B0 


2BQ  2B, 

— :L  \c  +  _ — 

B  2  ’  7  B2 


«  kx  ,  -kx  , 

+  c1Qe  +  cne  +  c12 


(3.95) 


where  c10  ,  ,  c ^  are  constants  and  k  =  \/-  B2  . 

The  six  constants  c7  ,  Cg  ,  c^  ,  Cjq  ,  c^  and  c^  are  determined  from  the 
boundary  conditions  (3.  69),  (3.  70),  (3. 72),  (3,  74),  and  (3.  77). 

From  (3. 69),  for  vq  (0)  =  0  we  get. 


c10  +  C11  +  c12  ~  0 


(3.96) 


From  (3.69)  for  v  '  (0)  =  0  we  get, 

B  2B  2Bo 

("Bj  ‘  ^")C7  +  C9  +  kc10  '  kcll  =  0 
From  (3.  70), 

B6  •  c8  +  B7k2  •  c,0  +  Byk2  •  cu  =  -P 


(3.97) 


(3.98) 
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where. 


2B7Bo 

B6  -  -T1  +  B8 


Qtt,3  (g-  +  Oft) 
7  3  (Qi'i2  +  b8  ) 


B8  “  (^f  +  Q1C1>C 


From  (3. 72), 


?  \ca  9  -ka 

B6ac?  +  B6cg  +  B?k  e  *  c1Q  +  B?k~e  K  •  cu  =  P 


(3, 99) 


From  (3.  74), 


q  ka  3  -ka 

B9c7  +  (B10k  +  Buk')eKa  "  (B1()k  +  BjjK  )e  Ka  cu  =  V 


{3, 100) 


where, 


Brt  =  J  B 


9  "8  (B1  +  1) 


t 

2 


^  3  rl  *  ,  ^lll  1 

_  +  — J— 


2(B1  +  1) 


B10 


Gt(Bj  +  1) 


Bi; 


"IT  <Qiti2tBi 


T^l3) 


From  (3.  77), 


Pt 

B12ac7  +  B12C8  +  B13k  e  c10  +  B13k~e  C11  =  M1  " 


2 


(3  101) 


where 


B  =  -  _ _ _ 

12  (Bj  t  1) 


Qttft 


Qj*!  4 

2  '  Bir) 


The  above  six  simultaneous  linear  equations  can  be  solved  to  obtain  the  six 
constants  which  when  substituted  in  (3, 95)  would  give  vq  (x),  u^  (x)  can  then  be 
obtained  from  (3,  93). 

3.  6  Condition  for  Uniform  Shear  Stress  in  an  Antisymmetric  Joint 

If,  in  a  joint,  the  two  adherends  are  identical  and  the  moment  M?  is  equal  to 
-Mj  ,  then  the  joint  is  antisymmetric  about  the  middle  of  the  adhesive  layer,  i.e.  about 
the  x-axis,  In  such  a  case,  it  is  possible  to  require  a  uniform  distribution  of  shear 
stress  along  the  x-axis,  and  derive  an  expression  for  the  ratio  (  V/P)  as  a  necessary 
and  sufficient  condition  for  ibis  requirement.  It  is  required  that. 


(t  )  n=  — 

'  xy  y  =  0  a 


G(u.  +  v  ’)  =  — 
M  o  a 


(3.102) 
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Substituting  (3. 102)  in  (3. 66), 


-  (f-  +  awiY  +  Q1C12(-UP  ’  0 


u"  -  2P 

U1  ■  TJ7 


(3. 103) 


where. 


=  +  <W  +  QA2 


(3. 103)  yields, 

U1  =  'WT  x2  +  blx  +  ^2  (3«  104) 

b^  and  b£  are  constants  of  integration. 

(3. 102)  yields. 


V  *  -3a  '  x2  +  ¥  +  b2>  <3-I05> 

From  (3.69),  since  v  '  (0)  =  0, 
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The  above  expressions  for  v  '  and  satisfy  the  differential  equation  (3.67)  identically. 
From  (3. 70), 


(-%-  +  QjtjtJbj  +  Q1t12b1  =  -P 


or 


bl  = 


-P 


<nr  +  QiV  +  Qi h2  > 


p 

17 


(3. 108) 


From  (3. 74), 


Pt  Q^i c 

a  i'T' 


4Q1t13(^rP>  -  v  -  o 


or 


—  +  (Q 1 2t  +  4Qltl3  )-4-  ’  v  =  0 

a  v  n 1  3X11  '  a1'/ 


or 


-1=  -L  + 

P  a 


(CVl  1  +  3  ^1*1  ^ 


a(  +  Qjtjt  +  Q^2  ) 


(3. 109) 


Boundary  conditions  (3.  72)  and  (3.  77)  are  identically  satisfied.  Two  special  cases  can 
be  considered. 


Case  (1)  : 


«  1 


The  ratio  (V/P)  then  reduces  to 


t_ 

a 


tl<t+  3  tp 

a  (t  +  tj) 


Case  (2)  :  —  <  <  1 
tl 


In  this  case  we  have 


(3.110) 


(3-  111) 


It  can  be  seen  that  as  the  ratios  -j—  and  increase,  the  required  value  of 

1  yl 

(V/P)  also  increases.  Eq.  (3.  Ill)  thus  provides  the  minimum  value  of  V/P  for  which  a 
uniform  distribution  of  shear  stress  can  be  obtained. 

A  similar  attempt  at  requiring  the  tearing  stress  cr^  to  be  uniform  along  the 
x-axis  results  in  the  trival  case  of  P  =  V  =  Mj  =  M2  =  0  indicating  that  the  tearing 
stress  distribution  cannot  be  completely  uniform  under  any  given  load  configuration. 

3.  7  Solution  for  the  General  Case 

In  the  case  where  the  adherends  are  not  identical  and  where  m  and  n  are 
integers  greater  than  1,  the  differential  equations  still  remain  linear  and  homogeneous 
but  their  number  is  larger  and  therefore  their  analytical  solution  becomes  more  complex. 
These  differential  equations  are  therefore,  solved  numerically  using  the  method  of 
linear  conversion.  The  details  of  this  method  and  the  computer  program  for  IBM  3o0 
are  given  in  Appendix  A. 
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4.  EXPERIMENTAL  INVESTIGATION 


4. 1  Introduction 

The  main  purpose  of  the  experimental  work  was  to  check  the  analytical  results 
and  establish  confidence  in  the  theory.  The  analysis  considers  the  variation  of  stresses 
through  the  thickness  of  the  adhesive  and  treats  the  problem  as  one  of  plane  strain. 
Photoelasticity  presents  the  most  suitable  means  of  experimentally  determining  stresses 
over  an  area  in  a  plane  problem  and  was  therefore  employed  in  this  case.  In  order  to 
check  all  the  stress  distributions  and  results  that  the  analytical  work  gives  would 
require  the  construction  and  testing  of  many  complex  specimens.  For  the  limited 
purpose  of  checking  the  validity  of  the  theory  however,  testing  of  two  specimens  of 
varying  overlap  lengths  was  thought  to  be  sufficient.  In  order  to  obtain  a  reasonably 
good  fringe  pattern  for  the  adhesive,  it  was  necessary  that  the  adhesive  layer  be  suf¬ 
ficiently  thick.  For  this  reason  the  joints  ere  made  of  a  laver  of  a  photoelastic 
plastic  simulating  the  adhesive,  bonded  between  two  identical  1”  thick  metal  adherends. 
The  joints  were  made  1"  wide  to  approximate  the  condition  of  plane  strain.. 

One  important  result  that  the  theory  predicts  is  the  effect  of  the  ratio  V/P 
on  the  distribution  of  stress.  This  effect  was  tested  experimentally  by  subjecting  the 
specimens  to  loads  at  different  angles  with  respect  to  the  plane  of  the  adhesive  layer. 

4. 2  Material  and  Specimen  Preparation 

Araldite  6020  epoxy  was  used  to  simulate  the  adhesive  layer.  This  epoxy  is 
available  in  soL;d  form  and  can  be  machined  to  the  required  shape.  Araldite  6020 
was  preferred  over  the  other  photoelastic  plastics  for  two  main  reasons.  Firstly, 
compared  to  the  other  plastics,  Araldite  6020  has  a  lower  fringe  value  which  means  a 
good  number  of  fringes  can  be  obtained  with  small  loads  thus  improving  the  accuracy 
of  the  experimental  results.  Secondly  and  more  importantly,  Araldite  6020  can  be 
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easily  bonded  to  aluminum  with  a  thin  layer  of  adhesive  of  the  same  material,  setting 
at  room  temperature.  Araldite  6020  adhesive  sets  at  room  temperature  when  used  with 
Ciba  hardener  951.  This  aids  in  keeping  the  initial  stresses  to  a  minimum.  In  fact 
with  a  little  care,  the  joints  can  be  made  practically  free  of  initial  stresses.  This 
fact  is  important  because  the  presence  of  high  initial  stresses  would  render  the  model 
useless  for  a  photoelastic  investigation. 

Adherends  were  made  of  1"  square  2024  T4  aluminum  alloy.  In  order  to  pre¬ 
pare  the  aluminum  surface  for  adhesive  application,  the  adherends  were  treated  in  a 
hot  acid  bath  at  150  - 160°  F  for  about  10  - 12  minutes.  The  hot  acid  bath  was  made  with 
66  ml.  of  sulphuric  acid  (H2S04),  315  ml.  of  water  and  13. 5  gms.  of  sodium  dichromate 
(Na2Cr9C>7). 

Two  pieces  of  Araldite  6020  epoxy  were  machined,  one  to  the  dimensions, 

?"  x  1"  x  2"  for  the  2"  overlap  model  and  another  to  the  dimensions,  x  1"  x  3", 
for  the  3"  overlap  model.  Each  piece  was  bonded  first  to  one  adherend  by  means  of 
an  adhesive  consisting  of  10  pbw  of  Ciba  6020  and  1  pbw  of  Ciba  hardener  951.  The 
joint  was  allowed  to  set  at  room  temperature  for  24  hours.  The  other  surface  of  the 
epoxy  was  bonded  to  the  second  adherend  in  the  similar  manner  and  allowed  to  set  24 
hours.  The  joint  was  tested  immediately  after  it  was  hardened,  because  the  epoxy 
starts  absorbing  mositure  from  the  air  and  high  initial  stresses  develop.  The  dimen¬ 
sions  of  the  two  specimens  are  shown  in  Fig.  6. 

4. 3  Test  Equipment  and  Procedure 

Fig.  7  shows  the  photoelastic  bench.  A  circular  polariscope  with  a  mercury 
green  monochromatic  light  source  was  used.  The  model  was  photographed  in  light 
background.  The  model  was  loaded  by  means  of  dead  weights  at  the  end  of  a  lever  arm 
with  a  mechanical  advantage  of  4:1.  The  direction  of  the  load  with  respect  to  the  plane 
of  the  adhesive  layer  '"as  varied  by  rotating  the  model  to  the  desired  angle.  This  was 


51 


achieved  by  means  of  two  aluminum  plates  fastened  to  the  model,  one  at  each  end  (Fig.  8), 
each  plate  provided  nine  loading  points  allowing  the  tangent  of  the  angle  (or  the  ratio  V/P) 
to  vary  from  zero  to  0.  8  in  steps  of  0. 1. 

At  first  the  model  was  photographed  with  no  load.  For  the  photoelastic  tests  to 
be  reliable  it  was  necessary  that  nothing  in  excess  of  a  half  order  fringe  be  present  at 
zero  load.  The  model  was  then  gradually  loaded  to  R  =  P  +  V  =  3C9  lbs.  The  move¬ 
ment  of  the  fringe  was  observed  and  the  model  was  photographed. 

In  order  to  obtain  the  stress  distribution  from  the  fringe  patterns  photographed, 
it  was  necessary  that  the  Araldite  epoxy  be  calibrated  for  its  fringe  value.  A  circular 
disc  of  the  material  was  loaded  in  compression  and  the  fringe  order  at  the  center  of  the 
disc  was  obtained  by  using  Tardy's  method  of  compensation  (16).  This  procedure  was 

repeated  with  different  loads  and  a  calibration  curve  was  obtained  (Fig.  9).  Theoret- 

4P 

ically  the  shear  stress  at  the  center  of  the  disc  is  given  by,  t  =  _ c_  where,  P  = 

7Td 

compressive  load,  and  d  =  diameter  of  the  disc.  The  fringe  value,  f,  of  the  material 

4P 

was  then  calculated  from  f  =  c  where  N  is  the  fringe  order  at  the  center  of  the 

?rdN 

disc. 

The  modulus  of  elasticity  for  Araldite  6020  (E  =  4.  5  x  10^  psi)  was  taken  from 
Leven  (17).  This  value  was  used  for  the  theoretical  stress  analysis  of  the  models.  The 
results  are  discussed  in  the  following  chapter. 
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5.  ANALYTICAL  AND  EXPERIMENTAL  RESULTS 


5. 1  General 

The  behavior  of  a  lap  joint  depends  on  a  number  of  factors  including  the  relative 
properties  of  the  two  adherends.  The  case  of  identical  adherends  is  of  special  interest 
and  will  be  discussed  in  detail.  The  performance  of  the  non-identical -adherend -joint 
will  be  compared  with  that  of  the  identical- adherend -joint.  The  case  where  the  joint 
is  subjected  to  pure  moment,  will  also  be  discussed. 

5. 2  Effect  of  m  and  n 

The  displacements  u,  and  v,  h^.e  been  expressed  in  Eqs.  (3. 1)  and  (3.2)  as 
power  series  in  y.  For  an  exact  elastic  solution  m  and  n  should  approach  infinity. 

It  should,  however,  be  possible  to  terminate  the  series  after  a  certain  number  of  terms 
and  arrive  at  a  solution  reasonable  close  to  the  elastic  solution. 

Fig.  10  shows  the  shear  stress  distribution  at  the  mid -plane  (y  =  0  plane)  of  the 
adhesive  layer  with  varying  values  of  m.  At  m  =  1,  the  displacement  u  is  a  linear 
function  of  y,  while  at  m  =  2  and  3,  the  displacement  u  is  expressed  respectively  as 
a  quadratic  and  a  cubic  in  y.  Qualitatively  it  is  seen  that  in  all  the  three  cases,  the 
maximum  stress  occurs  at  the  joint  edges,  decreasing  rapidly  to  a  minimum  in  the 
middle.  It  is  also  seen  that  as  m  increases,  there  is  a  decrease  in  the  maximum 
stress  concentration.  The  effect  of  m  is,  however,  much  more  clear  in  Fig.  11  where 
the  distribution  of  interfacial  shear  stress  is  shown.  For  m  =  2  and  3,  the  highest 
stress  is  seen  to  occur  at  the  leading  corner,  there  is  a  rapid  decrease  in  the  stress 
towards  the  middle  of  the  joint  and  then  another  smaller  peak  occurs  near  but  not  at 
the  trailing  corner.  The  stress  is  in  fact,  seen  to  decrease  at  the  trailing  corner.  This 
qualitative  distribution  is  not  seen  in  m  =  1  curve.  Theoretically  the  stress  at  the 
trailing  corner  should  be  expected  to  be  zero  whereas  the  singularity  at  the  leading 


.. ^ 
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comer  should  produce  a  high  shear  stress  at  that  comer.  It  is  interesting  to  note  that 
rising  values  of  m  produce  a  higher  stress  at  the  leading  comer.  This  indicates  that 
if  m  is  made  sufficiently  large  the  ideal  solution  may  be  approached.  The  real  mate¬ 
rials,  however,  do  not  behave  as  perfectly  elastic  materials  and  some  plastic  deforma¬ 
tion  may  occur  at  the  leading  comers,  the  extent  of  such  deformation  would  depend  on 
the  properties  of  the  adhesive  material.  The  stress  distribution  obtained  for  m  =  2  or 
3  can  be  expected  to  be  reasonably  close  to  the  actual  case  both  qualitatively  and 
quantitatively. 

The  distribution  of  a  at  the  mid -plane  and  at  the  interfacial  plane  is  shown  in 
Figs.  12  and  13.  At  the  interface  the  three  curves  are  in  good  agreement,  buc  at  the 
mid-plane,  m  =  2  and  3  show  cr  to  be  much  lower.  This  indicates  that  a  varies 

y  y 

appreciably  through  the  thickness  of  the  adhesive.  Similar  conclusions  can  be  drawn 
about  the  longitudinal  stress  crx  shown  in  Figs.  14  and  15. 

It  can  be  concluded  from  these  results  that  for  m  equal  to  2  or  higher  the 
qualitative  and  quantitative  distribution  of  stress  is  better,  compared  to  the  distri¬ 
bution  obtained  for  m  equal  to  1.  Also,  the  quantitative  results  for  m  equal  to  2, 
can  be  considered  sufficiently  close  to  the  elastic  solution. 

In  all  cases  discussed  above,  the  value  of  n  was  kept  constant  at  1.  The  thick¬ 
ness,  t,  of  the  adhesive  is  small  compared  to  the  overlap  length,  a.  For  this  reason 
the  numerical  integration  of  the  differential  equations  for  n  greater  than  1,  becomes 
extremely  inaccurate.  In  order  to  obtain  reasonably  good  results,  the  step  size  for 
integration  has  to  be  made  extremely  small.  The  computer  time  and  cost  involved 
would  in  such  a  case  be  high.  It  must,  however,  be  noted  that  for  the  very  fact  that 
t/a  is  small,  the  displacement  v  can  be  assumed  to  be  approximately  a  linear  function 
in  y,  and  n  can  therefore,  be  taken  as  1. 
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5. 3  Comparison  with  Experiment 

In  Fig,  16  are  shown  the  fringe  patterns  for  the  two  models  at  no  load.  It  is 
seen  that  the  zero  and  half  order  fringes  are  present.  Under  load,  fringes  of  up  to 
eighth  order  are  found.  The  models  can  be  considered  to  be  adequately  free  of  initial 
stresses. 

The  experimental  results  were  obtained  in  the  form  of  photoelastic  fringe  pat¬ 
terns.  The  computer  was  programmed  to  convert  the  maximum  shear  stress  at  each 
point  in  the  adhesive  layer,  into  the  corresponding  fringe  order.  The  comparison  of 
these  results  is  shown  in  Figs.  17  and  18  for  the  2  inch  specimen  and  in  Figs.  19  and  20 
for  the  3  inch  specimen.  The  theoretical  results  shown  in  these  figures  are  for  m  =  2 
and  n  =  1.  At  the  mid-plane  (Figs,  17  and  19),  the  theoretical  and  experimental  results 
are  in  good  agreement  except  at  the  two  ends,  where,  the  experimental  results  show  a 
rapid  decrease  in  stress.  At  the  interfacial  plane  (y  =  t/2),  again,  the  two  results  are 
in  excellent  agreement  except  at  the  trailing  end,  where  the  experimental  results  show 
lower  stress  concentrations.  It  can  be  argued  here  again,  as  has  been  done  in  dele 
5,2,  that  as  m  is  made  greater,  the  theoretical  stress  distribution  can  be  expected 
to  approach  the  experimental  pattern. 

The  analysis  predicts  that  the  ratio  V/P  has  a  significant  effect  on  the  stress 
distribution,  and  in  fact  for  a  certain  value  of  V/P  giver.  by  Eq.  (3. 109),  the  shear 
stress  distribution  can  be  made  uniform.  This  value  of  V/P  turns  out  to  be  0,  76  for 
the  2  inch  specimen  and  0.  51  for  the  3  inch  specimen.  Experimental  results  for 
V/P  =  0.  8  for  the  2  inch  specimen  are  compared  in  Figs.  21  and  22  with  the  corresponding 
theoretical  results.  At  the  mid-plane  (Fig,  21),  it  is  again  seen  that  the  two  results 
are  in  good  agreement  except  at  the  ends  where  the  experimental  stress  decreases 
rapidly.  At  the  interfacial  plane  (Fig.  22),  the  experimental  stress  shows  a  rapid 
decrease  at  the  trailing  edge  but  is  otherwise,  in  good  agreement  with  th  theory.  Simi- 


55 


lar  results  are  show  in  Figs.  23  and  24,  for  the  3  inch  specimen  for  V/P  =  0. 60.  It 
must  be  emphasized  here,  again,  that  the  theoretical  results  are  for  m  =  2,  and  can  be 
expected  to  approach  the  experimental  results  for  greater  values  of  m. 

In  order  to  compare  the  theoretical  and  experimental  results  over  the  whole  area 
of  the  adhesive  layer,  the  computer  was  programmed  to  plot  theoretical  fringe  patterns 
for  m  =  2.  These  fringe  patterns  are  compared  against  the  experimental  fringe  patterns 
in  Figs.  25  to  27  for  the  2  inch  specimen  and  Figs.  28  to  30  for  the  3  inch  specimen. 
Similar  conclusions  can  be  drawn  from  these  results  as  have  been  discussed  above. 

5. 4  Case  of  Identical  Adherends 

In  this  section,  the  analytical  results  for  a  joint  with  identical  adherends  and 
antisymmetric  loading  will  be  discussed.  Theoretically  it  has  been  found  that  the  stresses 
in  the  adhesive  layer  depend  on  the  factors,  E^/E  ,  tj/t  ,  tj/a  ,  V/P  and  P/a.  The 
stresses  can  be  expressed  in  dimensionless  form  by  dividing  them  by  the  mean  shear 
stress  P/a.  The  number  of  joint  parameters  then  reduces  to  four,  E^/E  ,  t^/t,  t^/a 
and  V/P, 

The  condition  for  a  uniform  distribution  of  shear  stress  at  the  mid-plane  of  the 
adhesive  layer  is  given  by  Eq.  (3. 109).  Let  Vf  be  a  factor  defined  as  follows 


Vf  = 


V 

P 


Q1t12(t  +  4- 1 


t  + 


3  T 


W 


(5.1) 


Then  the  condition  for  uniform  distribution  of  stress  is,  V^.  =  1.  The  factor  Vf  can 
therefore,  be  used  as  a  parameter  of  comparison.  The  joint  parameters  then  are, 
F^/E  ,  tj/t  ,  ^/a  ,  Vf. 

Let  and  be  parameters  defined  as  follows. 


S 


f 


(5.2) 
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(5.3) 


Tas  factor  S,  is  similar  to  /yfD  of  Vc lice r sen  (1)  ami  of  Goland  and  Reissner  (2). 

i  tj 

Tie  sfeear  stress  can  be  expressed  as  a  function  of  Sf. 


(t  j 

'  xvmax 


av* 


=  F,  (Sf ,  Vf,  Ej/E  ,  tj/t) 


(5.4) 


tafsere. 


av 


=  average  shear  stress  =  P/a 


(5,5) 


The  factor  N.  is  similar  to  the  parameter  X  of  Goland  and  Reissner  (2).  The  normal 


stresses  can  be  expressed  as  functions  of 


(o  ) 
vmax 


av 


=  F2(Nf,  Vf ,  Ej/E  ,  tj/t) 


(5.6) 


(c  ) 

X,m3X  *  -7—  =  F_  (N, ,  V  ,  E.  /E  ,  t./T) 
‘av  h  o  f  i  l  1 


(5,  7) 


In  the  following  discussion  on  the  shear  and  normal  stresses  in  the  adhesive 
layer,  the  theoretical  data  has  been  calculated  for  m  =  n  =  i  and  n  =  =  0.30,  It 

has  been  discussed  in  Article  5. 2,  that  m  =  2  gives  better  results  from  the  viewpoint 
of  accuracy  of  stresses,  m  =  1,  however,  provides  a  simpler  and  adequate  means  for 
comparison  of  joints, 

5.4.1  Shear  Stress:  Figs.  31  and  32  show  the  shear  stress  concentration  factor 
as  a  function  of  the  parameter  tj/t.  Curves  are  shown  for  various  values  of  E^/E  and 
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Sj  ,  keeping  constant  at  zero.  Figs.  33  and  34  show  the  same  curves  for  =  2. 
it  can  be  concluded  from  these  diagrams,  that  for  t^/t  2  10  and  E^/E  2  10,  the  shear 
stress  concentration  factor  is  independent  of  t^/t  and  E^/E,  and  is  a  function  of 
and  alone. 

From  design  point  of  view,  the  maximum  shearing  stress,  'rmax  »  is  even  more 
important  and  is  shown  in  Figs.  35  and  36  as  a  function  of  t^/t.  For  t^/t  2  10  and 
Ej/E  2  20,  the  maximum  shearing  stress  can  be  considered  as  a  function  of  and 
V,  only. 

From  the  above  discussion,  it  is  observed  that  the  shear  stress  concentration 
factor  and  the  maximum  shearing  stress  concentraion  factor  are  independent  of  tj  /t 
and  Ej/E  for  t^/t  2  10  and  E^/E  2  20,  and  are  functions  of  and  only.  It 
must  be  mentioned  here  that  the  peak  stress,  whether  positive  or  negative,  is  found  to 
occur  at  the  leading  comer  when  stresses  at  the  interfacial  plane  are  being  considered, 
and  at  the  two  ends  when  stresses  at  the  mid-plane  are  being  considered. 

In  Fig.  37  is  shown  the  shear  stress  concentration  factor  at  the  mid -plane  of 
the  adhesive  layer  as  a  function  of  and  V^.  It  is  interesting  to  note  that  the  shear 
stress  varies  approximately  linearly  with  S^,  and  also  with  V^.  For  >  6  the  shear 
stress  concentration  factor  can  be  approximately  expressed  as  follows. 


-  1  +  4<S(-1)(1-V  >. 

av 


at  y  =  0  plane  (5. 8) 


Fig.  38  shows  shear  stress  at  the  plane  of  the  interface,  as  a  function  of  Sf  and 
V^.  Again,  the  shear  stress  is  approximately  a  linear  function  of  and  .  The 
approximate  formula  in  this  case  is. 


(T  ) 
xy  max 
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av 
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i  +  -Z-  (Sf  - 1)  ( l  -  vf), 


at  v 
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(5,9) 
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Tmax  *s  s^ovm  as  a  Action  of  and  in  Fig.  39.  The  maximum  shearing 
stress  can  be  approximately  expressed  by  the  following  two  equations. 

For  ^  1  and  >  4; 

=  1  +  T(4I  '  V  <Sf  -  (5. iO) 

av  _ 

i 

For  Vj  -  1. 5  and  >'  5 
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)  (S.  -  1) 


(5.11) 


5.4.2  Normal  Stresses:  The  maximum  tearirig  stress,  (a  )  , 

-  &  \  y'max 

Figs.  40  and  41  as  a  function  of  the  parameters  t^/t  ,  E^/E  ,  and  „ 

J  1 


is  shown  in 
It  is  observed 


that  for  tj/t  ^  10  and  E^/E  £  10,  the  tearing  stress  factor  is  independent  of  the  para¬ 
meters  tj/t  and  E^/E  ,  and  is  a  function  of  Nf  and  Vf  alone. 

Figs.  42  and  43  show  that  the  variation  of  the  stress  (a  )  is  qpite  appreciable 

x  max 

for  Ej/E  <  20.  .However,  it  can  be  concluded  from  these  diagrams  that,  for  E^/E  -  20 
and  t./t  -  10,  (o)  is  a  function  of  the  parameters  Nf  and  V.  only. 

The  maximum  normal  stress  pt  any  given  point  acts  on  one  of  the  principal 
planes.  Fig.  44  shows  the  maximum  principal  stress,  omax  ,  as  a  function  of  the 
parameters  t^/t  ,  E^/E  ,  and  Vf  .  Again,  it  can  be  concluded  that,  for  tj/t  -  10, 
and  E./E  -20,  a  can  be  considered  to  be  a  function  of  two  parameters,  N,  and 
Vf  only. 

1  1 

The  maximum  tearing  stress  is  shown  in  Figs.  45  and  46  as  a  function  of  Nf 

and  Vf  .  These  diagrams  are  drawn  for  t^/t  =  10,  and  E^/E  =  20.  As  has  been 

discussed  atx  ve,  these  graphs  would  be  applicable  for  all  joints  where  t^/t  -  10  and 

E^/E  -  20,  It  must  be  noted  here  that  this  range  is  the  one  that  is  most  frequently 
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encountered  in  practice.  The  data  shown  in  Figs.  45  and  46  are  therefore,  applicable 

to  the  majority  of  the  practical  lap  joints. 

Similar  graphs  are  shown  in  Figs.  47  and  48  for  ( c  )  ,  and  in  Fig.  49  for 

x  max 

cmax*  ^  *S  °bservec*  t^at  f°r  a  given  V^,  the  normal  stress  factors  vary  approximately 

as  square  functions  of  N^.  This  is  more  clearly  seen  in  Figs.  50  through  54,  where 

2 

the  stress  factors  are  shown  as  functions  of  .  These  graphs  consist  of  almost 
straight  lines.  The  stresses  can  '  j  approximately  expressed  as  follows. 
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Y~  =  (0.  85  -  0. 95  Vf)  Nj  ,  at  y  =  +  — 
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Now,  Eqs.  (5.  2)  and  (5. 3)  yield, 


N  2  =  - a  •  S 
f 


Substituting  (5, 17)  in  (5. 12)  through  (5. 15), 


l.y?.max  =  (1.325  -  1,6  V)S 
av  ' 


at  y  =  C 


(5.17) 


(5.18) 
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max  =  (2.0  -  2.2  V)  Sf, 
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at  y  =  +  -2 


(5.22) 


Thus  it  is  seen  that  the  normal  stresses  are  approximately  linear  functions  of 
Vj  and  Sf  .  Eqs.  (5. 18)  through  (5. 22)  show  another  interesting  development  if  both 
sides  of  these  equations  are  multiplied  by  t^/a  ,  which  yields. 


(°y^max  *1  _  „  oc  .  a  m  \  .  ^  A 
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Since  is  a  linear  function  of  the  overlap  length,  a  (Eq,  5,1),  it  can  be  con¬ 
cluded  that  the  normal  stresses  are  also  linear  functions  of  the  overlap  length. 

5.4.3  Effect  of  :  It  has  been  observed  that  is  an  important  factor  con¬ 
trolling  the  stress  concentrations.  It  would  therefore,  be  of  interest  to  look  at  the  stress 
distributions  along  the  overlap  length  for  different  values  of  .  Fig,  55  shows  the 
shear  stress  distribution  at  the  interfacial  plane  for  three  values  of  V{.  The  maximum 
shear  stress  occurs  at  the  leading  comer  and  is  in  the  positive  direction  for  =  0; 
and  in  the  negative  direction  for  =  2.  At  =  1,  the  shear  stress  at  the  inter  - 
facial  plane  is  nearly  uniform.  Figs.  56  and  57  show  the  distribution  of  and  ox  . 
Again  for  =  0  and  2,  the  maximum  stress  occurs  at  the  leading  comer  and  is 
positive  for  =  0,  and  negative  for  =  2  .  At  =  1,  the  normal  stresses  are 
nearly  uniform  over  much  of  the  overlap  length  except  near  the  leading  corner,  where  a 
negative  maximum  occurs, 

5.4.4  Design  of  Lap  Joint:  The  importance  of  the  factors  Sf  and  Vf  in 
relation  to  the  stresses  in  the  adhesive  layer,  has  been  demonstrated.  The  stresses 
are  greatly  reduced  as  the  factor  is  brought  closer  to  1.  This  fact  can  be  usefully 
employed  in  designing  lap  joints.  In  practice,  lap  joints  are  used  to  fasten  two  sheets 
subjected  to  pull.  The  pull  can  be  applied  at  any  angle  0  along  the  length  of  the  adherends 
(Fig.  58a),  The  equilibrium  of  the  joint  requires  that  the  pull  at  the  two  ends  be  colineai. 
In  this  case  the  ratio  V/P  is  given  by, 

Y  (t,  +  t) 

_V_  =  tan  0  =  —1-r-  (5.28) 

&  is  the  length  defined  in  Fig.  58a,  The  minimum  value  of  !>  that  can  be  used,  is 
equal  to  the  overlap  length,  a.  The  maximum  value  of  the  ratio  V/P  therefore,  in 
this  case  is  given  by, 
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V  _  (tI  +  *> 

P  a 

Further,  if  t/t^  is  small  then  the  maximum  value  is. 


P  a 


(5.29) 


(5.30) 


Eq.  (3.  Ill)  gives  the  required  value  of  the  ratio  V/P  for  uniform  distribution  of  stress 
as, 

V  _  4  .  *1 

P  3  a 


Thus  the  maximum  value  of  V/P  obtainable  by  means  of  simple  bending  of  adherends. 
is  less  than  that  required  for  uniform  shear  stress  distribution.  In  other  words,  the 
factor  Vf  for  the  joint  shown  in  Fig.  58a,  is  less  than  1.  The  value  of  this  factor 

is  approximately  0.  75  for  ^  =  a.  This  value  is  fairly  close  to  1,  and  the  stress  curves 
show  that  this  design  provides  an  excellent  improvement  over  the  conventional  straight 
adherend  design.  It  is,  however,  possible  to  further  improve  the  design  by  applying 
the  pull  at  angle  to  the  adherend  (Fig.  58b)  and  requiring  that  the  condition  for  =  1 
be  satisfied. 

5.  5  Case  of  Pure  Moment 

A  lap  joint  is  under  pure  moment  when  P  =  0,  V  =  0,  and  .  In  Figs.  59, 

60,  and  61  are  shown  the  stresses  in  a  joint  subjected  to  pure  moment.  The  stresses 
were  determined  for  m  =  2.  The  peak  stresses  are  found  to  occur  at  the  two  ends  of 
the  joint.  The  highest  stress  occurs  at  the  leading  corner,  as  would  indeed  be  expected. 
The  stresses  are  positive  at  that  end  of  the  joint  where  the  moment  produces  tensile 
stresses  in  the  adherend  at  the  interfacial  plane.  The  peak  stresses  are  negative  at  the 


other  end. 


5. 6  Non-Identical  Adherends 


The  case  of  identical  adherends  has  been  discussed  in  detail.  The  performance 
of  the  non-identical  adherend  joint  will  now,  be  compared  with  that  of  the  identical 
adherend  joint.  The  two  adherends  in  a  joint,  can  differ  in  their  material  properties 
as  well  as  in  their  thicknesses.  Theoretical  results  were  obtained  for  a  case  where 
t2  and  E2  were  varied  while  holding  t^  and  constant.  The  case  of  t2  =  t^  and 
E2  =  Ej  provided  the  identical  adherend  joint  as  the  base  for  comparison.  The  ratios 
t2  /tj  and  E2/Ej  were  increased  from  their  initial  value  of  1.  The  stress  at  the 
leading  comer  of  adherend  (I),  which  is  also  the  highest  stress  in  the  joint,  was  obtained,* 
This  stress  is  presented  in  Figs.  62  through  67,  as  a  multiple  of  the  stress  for  the  base 
joint.  The  results  are  shown  for  two  different  loading  conditions,  (i)  V  =  0,  and 
(ii)  Mj  =  M2  =  0.  It  was  found,  in  both  cases,  that  for  E^E^  =  1,  10,  and  100,  the 
results  were  almost  identical  indicating  that  the  factor  E^E^  has  no  appreciable  effect 
on  the  maximum  stress  in  the  joint.  The  factor  t? /t^  however,  has  a  marked  effect  on 
the  maximum  stress,  which  increases  as  t2/t^  increases.  The  ratio  of  this  increase 
is  much  greater  in  the  case  where  V  =  0,  as  compared  to  the  case  where  =  M2  =  0, 
The  moments  and  M2  when  V  =  0,  are  given  by. 


M 


1 


(5, 30) 


Thus  Mj  is  proportional  to  t9  ,  which  explains  the  rapid  increase  in  stresses  as  t2 
increases.  This  also  indicates  that  the  increase  in  M^,  rather  than  the  increase  in 
stiffness  of  the  adherend  (II),  is  responsible  for  the  increase  in  stresses.  One  important 


*  For  t9/tj  -  1  ,  and  E2/E^  ^  1  ,  the  maximum  stress  would  occur  at  the  adherend  (II) 
interface  at  its  leading  comer. 


conclusion  from  these  results,  is  that  the  joint  with  identical  adherends  is  the  most 
efficient  joint  and  that  any  difference  in  the  thicknesses  of  the  two  joints  teiJs  to  increase 


the  maximum  stress. 
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6„  SUMMARY  AND  CONCLUSIONS 

An  elastic  analysis  of  the  lap  joint  has  been  presented.  The  two  adherer.ds  are 
allowed  to  differ  in  their  thicknesses  and  mechanical  properties.  The  joint  is  subjected 
to  a  general  loading,  consisting  of  tension,  shear  force  and  bending  moments.  The 
stresses  in  the  adhesive  layer  are  determined  treating  the  problem  as  one  of  plane 
strain.  The  analysis  is  therefore,  applicable  to  joints  with  large  widths.  The  main 
limitation  of  the  theory  is  that  it  requires  the  adhesive  material  to  be  isotropic  and 
linearly  elastic. 

In  order  to  check  the  analytical  results,  and  to  establish  confidence  in  the  theory, 
two  specimens  of  lap  joint  of  varying  overlap  lengths  were  made  and  tested  photoelas- 
tically.  Excellent  agreement  was  found  between  the  theoretical  and  the  experimental 
results,  thus  establishing  validity  of  the  theory. 

The  case  of  a  joint  with  identical  adherends  and  subjected  to  antisymmetric 
loading,  was  analytically  studied  in  detail.  The  following  conclusions  can  be  drawn 
for  this  joint. 

(1)  The  maximum  shear  stress  and  the  maximum  normal  stresses  are 
approximately  linear  functions  of  two  parameters,  and  .  Simplified  formulas 
giving  the  maximum  stresses  as  functions  of  and  Vf  are  obtained. 

(2)  The  maximum  normal  stresses  are  linearly  dependent  on  the  overlap 

length. 

(3)  A  nearly  uniform  distribution  of  the  shear  stress  can  be  obtained  for 
Vj  =  1.  The  normal  stresses,  although  not  uniformly  distributed,  are  minimized  in 
this  case.  Thus  =  1  provides  the  most  efficient  joint. 

(4)  Design  of  the  lap  joint  based  on  the  above  results,  has  been  discussed  and 


recommendations  made. 
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The  analytical  results  for  a  joint  with  identical  adherends,  and  subjected  to  pure 
moment,  show  high  stress  concentrations  at  the  two  leading  comers.  The  adherend  that 
is  subjected  to  an  anticlockwise  bending  moment  ha,*,  at  its  leading  comer,  high  positive 
stresses.  The  stresses  at  the  other  leading  comer  are  negative. 

Following  conclusions  can  be  drawn  from  the  analysis  of  joints  with  non-identical 
adherends, 

(1)  The  difference  between  the  modulii  of  elasticity  of  the  two  adherends,  has 
no  appreciable  effect  on  the  maximum  stress  in  the  joint. 

(2)  An  increase  in  the  thickness  of  one  adherend  as  compared  to  that  of  the 
other,  results  in  an  increase  in  the  maximum  stress. 

(3)  Considering  variation  in  the  adherends,  the  joint  with  identical  adherends 
is  the  most  efficient. 
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APPENDIX  A 

SOLUTION  OF  DIFFERENTIAL  EQUATIONS 
IN  THE  GENERAL  CASE 

Equations  (3. 51)  through  (3. 54)  are  linear,  homogeneous,  ordinary  differential 
equations.  The  boundary  conditions  are  partly  at  one  end  and  partly  at  the  other.  The 
problem  can  be  transformed  into  an  initial  value  problem  by  using  linear  conversion. 
The  equations  are,  at  first,  reduced  to  first  order  differential  equations  by  the  fol¬ 
lowing  substitution. 


where, 

i  0,  1,  2,  m 

j  -  0,  1,  2,  n 

ij  =  i  +  (m  +  1) 

i2  =  i  +  2  (m  +  1) 


j3  =  j  +  3  (m  +  1) 


TO 


i4  =  j  •=■  3(123-5-  S)  -5-  (ail) 

j-  =  j  -5-  3(na-S-l)  -5-  2(a-5-l) 

L  =  j  -5-  3(121-5-1)  -5-  3(e -3- 1) 


Eteffluas  (3.51)  ami  (3. 52)  reduce  »  (sa  -r  1)  algebraic  equations,  given  by. 


m 

T" 


i  =  0 


EdB 

,  k)  w  *  ^ 

1  i=0 


b12  t5'  k)  % 


-  =  bp(i.  k) 

j  =  0 


a 

■v  2  Ijl  .  (j,  k)  w.  =0  for  k  =  0,  1,  2, . . . .  m  (A  2) 

3=0  iS>  J6 

is^liere, 

Ugj  {«,  C)  —  0,  i  =0y  I,  2,  ........  m 

0^(0,  k)  =  0,  k  =  0,  1,  2, . .  m 

b.&  (i,  t)  =  -Gik  (  J-)1^  k  '  1  o.k  ,  i  =  1,  2 . .  m  k  =  !,  2 . m 

by,  <«.  0)  =  (4-)1  '  %^io  +  Qj1!  -  (  l)1  1  =  0,  1.  2 . .  m 

w 

b12  o*  »  *  -<T>’  *  4  ;  -r  “it  ’  •Vi  +  Q2t2(  l,i  * *  j 

4» 

i  *■  0,  if  2f '  •  *  *  f  in »  ^  2,  • .  •  •  t  Hi 

bi4  <0,  0)  -  0 


bJ4  <J»  *?>  =  '  l3 


;o 


1  =  1- 


71 


b14  (0,  k)  =  -Gk  (\  f  ^Qk 


k  =  1  2 

•  *  »  <i 


m 


bl4(j,  k)  =  (\j  -Gk)(-|-)  j  +  k  ^jk,  j  =  1,  2,  .....  n 

y  =  1,  2,  .  m 


b160,  0)  =  f  -Q^2  -r  (-1^  Q2t22]  , 

bl60,k)  =  [QiV  *  H)j"k  %h 


j  0,  If  2,  • « « « i  n 

^  j  0^  If  2,  » f ;  j  •  n 

k  =  !.  2t  m 


Equations  (3. 53)  and  (3. 54)  are  reduced  to  (n  -f  1)  first  order  differential 
equations,  given  by. 


in  n  rn 

£  b  9(i,  k)w  '  +  £  b^,(j.  k)w  '  =  s  b„  0-  k)w 

i  =  0  l2  j  =  0  J6  i  =  0  *1 


n 


n 


H-  S  b  o.  S)«'  4-  2  b,5(j.  k)«j_ 

j  =  0  J3  j  =  0 


=  0,  k  =  0,  1,  2 . .  n 

(A.  3) 


vhere. 


b2,  (0,  0) 


=  C 


b21  (i.  0)  =  Q  (4-)1  A 


10 


i  -  1,  2, . m 


b?.  (i,  k)  =  (Gi  *  Xk)(4>1  +  k  A 


ik 


i  -  0.  1,  2. 

k  -  1,  2. 


m 


(l,  * 


«-W 


'Vf 


<  -!»’  Q2T2~ j 


i  •  •  >  I  2 
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_  j  ,  t  J  +  k 


b22  (i,  k)  =  {  (^-) 


-QX2  +  (-l)i  +  k  Q0t,2 


^2  2 


i  =  0,  1,  2, 
k  =  1,  2,  . . 


. . . ,  m 

i«i«y  n 


(0,  k)  =  0,  k  =  0,  1,  2,  n 


b23  Cj»  °)  —  0,  j  =  0,  1.  2, - -  n 


b23(j,  k)  =  -Qjk(-§-)j  +  k  _1  6jk,  j  =  1,  2,  ,...,n 


k  =  1,  2, 


n 


b25a,o)  =  -G(-i-)j+1ojo, 


j  Of  If  2,  •  •  •  •  f  n 


b25a  k)  =  G(~)j'hk+l  ,ajk  j  =  0,  1,  2,  .....  n;  k  =  1,  2,  .....  n 


b26  0.  °)  =  -  -J  <T  9 


Qft3  +  (-l)j  Q2t23 


f  j  Of  If  2,  « •  •  *  f  n 


b26  <i*  k>  =  -y(T) 


=  i  _L  ^  Int^ 


i-rk 


Qjtf  +  (-!)•  Q2t2 


j  0>  If  2,  •  •  •  •  t  n 
k  =  1,  2 . .  n 


Equations  (A.  2)  can  be  differentiated  once  to  give  additional  (m  -f  1)  first  order 
differential  equations  as  follows. 


m 


m 


£  b  (i,  k> »  '  +  _  Z  b  0.  k)  w.  2  bio  <j-  k>  "•< 

x  =  0  2  j  -  0  Jo  i  =  0  1 


-  b.  4  (j»  k)  w 

j  =  0  J5 


for  k  =  0,  1,  2 . .  ir 


(A.  4) 


Equations  (A.  3)  and  (A.  4)  can  be  treated  as  st^vltaneous  equations  in  and 


w.  *  and  can  be  solved  as  such,  to  obtain  the  following  system  of  equations. 

J6 

wi,'*Vwo'wi; . >'V  </L5) 


and 

w  ’  =  f  (w ,  w  , . .  w  )  (A.  6) 

J6  J6  u  q 

where,  q  =  3m  +  4n  +  6 

Also,  equations  (A.l)  can  be  written  in  the  form, 


'V  =  wi  =  fi<'V  wr  •• 


•  V 


v’i  ’  =  W  =  f  (w  ,  w„ 
i.  i2  1^  0  1 


’  '  q) 


V/.  ’  -  w.  =  f.  (w_,  w. 

J4  13  0  ■ 


"  'V  > 


(A.  7) 


\  =  "j5  =  fi4  (wo’  W1 . 'V 


w:  ’  =  w.  =  f.  (wn,  w. 


j5  J6  j5  '”0'  T 


,  v:  ) 

q 


Now,  the  method  of  linear  conversion  can  be  applied  to  obtain  w.(0),  for 
i  =  0,  1,  2, .  q.  The  following  equations  are  used. 


q 

w  (a)  =  -  a  w  (0)  4-  p  1  =  0,  1,  2, . .  q  <A.o> 

1  j  =  ft  }  1 


In  order  to  determine  the  con  starts  ar.d  the  following  steps  are  fol¬ 


lowed. 


1 


s 

I 

I 

I 


1 
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(1)  Set  wq(0)  =  w^(0)  =  w2(0)  =  .  =  w^(O)  =  0  and  integrate  equations 

(A.  5),  (A.  6),  and  (A.  7)  using  Runge-Kutta  fourth  order  formula.  This  gives, 

=  w.(a),  i  =  0,  1,  2, - ...»  q 

In  this  case  all  the  pi  are  zero. 

(2)  Set  w  (0)  =  1,  w.(0)  =  w9(0)  = .  =  w  (0)  =  0  and  integrate 

equations  (A.  5),  (A.  6),  and  (A.  7).  This  yields, 

a  j  =  w.(a)  -  P£  =  w.(a),  i  =  0,  1,  2, . .  q. 

This  step  is  repeated  for  w^(0)  =  1  and  all  other  zero,  etc. ,  until  all  a.  ^  are 
determined. 

(3)  (4n  +  2m  +  6)  boundary  conditions  are  known.  Equations  (A.  2)  supply 
another  (m  +  1)  boundary  conditions,  so  that  a  total  of  (q  +  1)  boundary  conditions  are 
known.  Equations  (A.  8)  provide  an  additional  set  of  (q  +  1)  equations.  These  (2q  +  2) 
equations  can  be  solved  for  (2q  +  2)  initial  and  final  conditions. 

(4)  Using  the  initial  conditions  from  step  3,  integrate  equations  (A.  5),  (A.  6), 
and  (A.  7)  a  final  time  to  get  the  answers.  As  a  check,  w.(a)  obtained  via  this  inte¬ 
gration  must  agree  with  those  obtained  previously  in  step  3. 

The  Runge-Kutta  fourth  order  formula  to  be  used  for  integration,  can  be  stated 
as  follow's. 

w*.  (x-rh)  =  w.  (x)  +  -l-  (a.  +  2b.  -f  2c.  +  d-  ) 
i  '  i  6  •  i  i  i  i' 


h  =  step  size 


r 


TV 


q 


where 
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a  a.  a 

bi  =  hf.  wQ  (x)  +  -y-,  w.  (x)  +  *2-  w  (x)  + 

b  b.  b 

ci  =  hfi  w0  (x)  +  ~T  »  wi  <x>  +  —  »  ••••••..  'vq  (x)  +  -|- 


di  =  hfi  Wo  (X)  +  co’  wi<x>  +  ci  »  •••••••»  wq(x)  +  cq 

The  computer  program  using  the  above  method  for  solving  the  differential 
equations,  is  given  in  Appendix  B. 


o  o  o  o  n  o 
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APPENDIX  B 

FORTRAN  PROGRAM  BASED  ON  IBM  SYSTEM  360/75 
FOR  THE  GENERAL  CASE 


Notation: 

E»E1 * E2  = 

E,  Ej  ,  E2  respectively 

D,D1,D2  = 

G,  Gj  ,  G2  respectively 

B,B1,B2  = 

Q/2  ,  Qj/2  ,  Q2  /2  respectively 

T,T1,T2  = 

t,  t^  ,  t2  respectively 

A  = 

Overlap  length,  a 

P»  V»  Cl *C2  * 

P,  V,  M.  ,  M2  respectively 

M  = 

m  + 1 

N  = 

n  +  1 

MM  = 

Number  of  data  points  along  x-axis,  (M  M-l) .  Ax  =  a 

NN  = 

Number  of  data  points  along  y-axis,  (NN-1).  Ay  =  t 

L= 

Integer  denoting  number  of  divisions  of  Ax  for  integration 

H  = 

Step  size  for  integration 

L*H*(KM-1.0) 

-  A,  i. e.  L '  H  =  Ax 

SIGMAX  = 

0 

X 

SIGMAY  = 

o 

y 

TAUXY  = 

T 

XV 

TAUHAX  = 

T 

max 

THE  PURPOSE 

OF  THIS  PROGRAM  IS  TO  DETERMINE  STRESSES  IN  T 

ADHESIVE  LAYER  OF  A  LAP  JOINT  OF  GIVEN  GEOMETRY  AND  LOADING. 


INPUT  =  EfEIfE2»0»DI*02»BfBl»B2»T*Tl»T2»PfV»Cl»AfHfH,N,HHfNN,l 
OUTPUT  =  SIGMAX,  SIGHAY,  TAUXY,  TAUHAX 
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IMPLICIT  REAL*8(A-H.,0-ZI 

DIMENSION  F(NW2,NW2) ,G( NW2 ) »B10 (M ,M ) , B12 ( M,H) , AA( MN,MN) , W(NW, MM ) 
l,B14IN,Mlt616(N,H),B22IM,N),B26NN,N),B2i(M,N) ,B23(N,N) ,B25( N,NS , 
2AU2(M,M) , AVl(NtM) ,AV3(N,M) ,CVl (N, N) ,CUO (M , N ) , DU1(M,N)»DV2(N»N) 
3AX( NW ) ,  BXINW) ,CX( NH ) ,DX (NW ) ,GX( NW2 , 1 ) 

4GGI MN» MN2 J ♦ ALFA (MN*  MNI *  BETA(MN*MN ) ,DELT ( MN, MN) »  SI GMAX(MM*NN ) 

5, SIGMAY I  MM, NN ) , TAUXYJ MM  fNN) , TAUMAXC  MM  »NN) 

MN=M+N 

MN2=M+2*N 

NW=3*M+4*N 

NW2=2*NW 

C2=Cl+V*A-P*(Tl/2.0+T2/2.0+T) 

BB=2.0*(B-D) 

C  DETERMINE  THE  COEFFICIENTS  FOR  THE  BOUNDARY  CONDITIONS  AND 
C  THE  DIFFERENTIAL  EQUATIONS 
821(1,1  )=0.0 
00  6  1=1, M 
DO  2  J=1,M 

AU2(I  ,J)=2.0*(  T/2.0)**(  I+J-2  )*(B+T/ 2.0*  Cl. 0+( -1.0) ++U+J-2)  )/ 

1< I+J-1J0)+B1*T1+B2*T2*(-1.0)**( I+J-2) ) 

2  CONTINUE 
DO  3  J=1,N 

AV3\ J»I )=( T/2 .0 )**( I+J-2 ) *(-Bl*Tl**2+( - 1. 0 ) ♦*( I +J-2 )*B2*T2**2) 
DU1( I,J)=(T/2.0)**( I+J-2l*(-Bl*Tl**2+(-1.0)**( I+J-2 )*B2*T2**2) 

3  CONTINUE 
AV1( 1,1 )=0.0 

B22 ( 1 , 1 )=-DUl (1,1) 

DO  4  J=2,N 

BETA!  I,  J)=(1.0-(-l.  0)**(  I+J-2 ))/(  I+J-2. 0) 

AV1U.I  )=BB  *(  J— 1 .0)  *(T/2. 0)**(  I+J-2  )*BETA(  I,J) 

B22( I » J)=0U1( I, J) 

4  CONTINUE 
610(1,1) =0. 0 
B10( 1 » I  1*0.0 
B12( T , 1 i=AU2( 1,1) 

DO  5  J=2»M 

B 1 2  (  I  ,'J  )=-AU2  ( I  ,  J  ) 

5  CONTINUE 

6  CONTINUE 

DO  10  I =2 , M 
DO  7  J=»2,M 

DELT( I»J)=(1.0-(-1.0)**< I+J-3 ) ) /( I+J-3.0) 

B10C I , J  )=— D=M I— 1.0) *( J-1.0)*( T/2. 0)**( I+J-3 )*DELT ( I » J) 

7  CONTINUE 

81411,1 )=-D*( T/2.0)**(I— 1)*(1 .0— {— 1.0 )**( I— 1) ) 

DO  8  J=1,N 

BETA{ I , J)=( 1 . O-C-l .0)** (I +J— 2 ) ) /( I+J-2. 0) 

CUO(I,  J)=D*(  I-1.0)*(T/2.0)**(  I+J-2)*BETA(I,J) 

B16( J, I )*-AV3( J,I ) 

8  CONTINUE 

B2 1 ( 1 , 1 )=— CUO  C I , 1 ) 

DO  9  J=2,N 

821 ( I »  J  )*CUO( I,J)-AV1(J,I) 

821(1, J)=-AV1(J,1) 
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jjf 

5 


I 


B14( J, I )=-B21 ( I  ,  J  ) 

9  CONTINUE 

10  CONTINUE 

DO  13  1=1, N 

cuod,  n=o.o 

B16(I *1)=AV3( 1,1) 

823(1, I)=0.0 
B23(I,U=0.0 
814(1, 1>=-AV1(I,1) 

DO  11  J=1,N 

ALFAU,  J|  =  (1.0-(-1.0)**(I+J-l)J/(  1+J-l.O) 

CV1( 1 ,Jl=(T/2»0)**( I*J-1)*D*ALFA( I,JI 

DV2( I , J) =2. 0*(T/2«0)**( l+J-2)*( B1 *T 1**3 /3«0+B2*T2**3*(- 1.0) 
1**( I+J-2) /3.0) 

11  CONTINUE 
826(1*1 )=-DV2( 1,1  ) 

825(1,1 )=-CVl (1,1) 

DO  12  J=2,N 
B26(I,J)=DV2(I,J) 

B25( I , JI=CVi( I, J) 

12  CONTINUE 

13  CONTINUE 

DO  15  1=2, N 
DO  14  J=2,N 

DELT( I,J)=(l,0-(-1.0J**(I*J-BI)/( I+J-3.0) 

B23( I,J)=-2.0*B*( I-1.0)*( J-1.0>*( T/2.0>**( I *J-3 »*DELT ( I ,J» 

14  CONTINUE 

15  CONTINUE 

C  SOLVE  SIMULTANEOUS  EQUATIONS  ( A. 31  AND  (A. 4)  OF  APPENDIX  A 

C  AA  ARE  THE  COEFFICIENTS  ON  THE  LKS  AND  GG  ON  THE  RHS 

DO  29  1=1, M 
DO  27  J=1 ? M 
AA(J,IJ=B12(I,J) 

GG(J,IJ=B10(I,J) 

27  CONTINUE 
DO  28  J=1,N 
J1=J+M 
J2~J1*N 

AA( J1,I  )=B22(I,J) 

GG(J1,I )=B21(I,J> 

AA(I,J1)=B16(J,I) 

GG(  I  j«Jl  1=0»  0 
GG( I , J2 )=814( J, I ) 

20  CONTINUE 

29  CONTINUE 

DO  31  1=1, N 
DO  30  J»1,N 
I 1=I*H 
Jl= J+M 
I2=I1*N 

AA(J1,I1»=B26(I,JJ 
GG( J1,I l)=823( I , J) 

GG(J1»I2)=B25( I , J I 

30  CONTINUE 

31  CONTINUE 


& 


o  o  o 
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CALL  SIMEQ  (AA,GG,MN,MN2, KOI 
IF(KO.EQ.l)  GO  TO  3000 
NOW  USE  LINEAR  CONVERSION 

F  AND  G  ARE  THE  COEFFICIENTS  ON  LHS  AND  RHS  RESPECTIVELY 
FIRST  SET  ALL  F  AND  G  TO  ZERO 
DO  26  J=i,NW2 
DO  25  1=1, NW2 
F( I ,J)=0.0 

25  CONTINUE 
G( J)=0.0 

26  CONTINUE 

C  NOW  DETERMINE  F  AND  G  THAT  ARE  NOT  ZERO 
DO  50  J=l , NW 
DO  35  1=1, NW 
W(I ,11=0.0 
35  CONTINUE 
W(J, 11=1.0 

C  INTEGRATE  EQUATIONS  ( A. 5 ) , < A. 6} , ( A .71 

CALL  INTEG  (GG,W,H, M, N, MM, NW, AX ,BX,CX,DX,  L ) 

DO  40  1=1, NW 
F ( I, J)=W< I, MM) 

40  CONTINUE 
JW= J+NW 
F(J,JWl=— 1.0 
50  CONTINUE 

C  USE  EQUATIONS  (A. 2) 

DO  80  J=1,M 
JW=J+NW 
DO  60  1=1, H 
I2=I+2*M 

F ( JW, l )=— B10( I  ,  J) 

F( JW, 121=812(1, J) 

60  CONTINUE 
DO  70  1=1, N 
14=1  -»-3*M+N 
I6=I4+2*N 

F(JW,I4J=-B14II,JI 
F( JW,I6)=B16tI,JJ 
70  CONTINUE 
80  CONTINUE 

C  NOW  USE  THE  BOUNDARY  CONDITIONS 
DO  400  1=1, M 
I1*I+H 
I2=I1+M 
I7=I+NW 
18=1 1*NW 
I9=I2+NW 
DO  90  J=2,H 
JW1=J«-NW+H 
Ft  JWl , II )=AU2C I , J  J 
GUwn=P*(-T/2.0)**IJ-l  J 
90  CONTINUE 

DO  100  J=1 , M 
JW2=J*NW+2*H 
FC JW2,18i=AU2(I,J) 


80 


GlJW2)=P*tT/2.9)**tJ-ll 
100  CONTINUE 

DO  200  J=2»N 

JW3=J+NW+3*N 

JW4=JW3+N 

Ft JW3f I )=CU0( ItJ) 

FI  JH3f  I2)=-DU1U,  J) 

GCJW3)  =  V*(- T/2. 0) ** ( J-l ) 

F ( JW4«  I 1 )=DU1 ( I ,J ) 

Gt JW4)=(-T/2.Q)**< J-l }*{ P*T2/2.0+C2) 
200  CONTINUE 

DO  300  J=1,N 

JW5=J+NW+2*N+3*M 

JW6=JW5+N 

Ft JW5,I7)=CU0tI,J> 

Ft JW5,I9)=-DUl<I,J) 
G(JW5)=V*tT/2.0>**<  J-l) 

Ft  JW6* I8)=DU1 1 1 » J } 

G( JW6)=-t T/2. 0)** (J-l )*tP*T 1/2. 0-^1 i 
3 00  CONTINUE 
400  CONTINUE 

DO  900  1=1, N 

I3=I+3*M 

I4=I3+N 

I5=I4+N 

16=1 5+N 

110=1 3+NW 

I11=I4+NW 

I12=I5+NW 

113=16+ NW 

DO  500  J=2,H 

JW 1= J+NW+H 

F ( JW1 » 1 3 ) =AV1 1 I , J ) 

Ft  JW1 » I5)=A V3 (I , J) 

500  CONTINUE 

DO  600  J=1,M 
JW2=  J+NW+2*M 
Ft JW2,I10)=AVltI, J) 

Ft  JW2,U2)  =  AV3(I,J) 

600  CONTINUE 

DO  700  J=2,N 
JH3=J+NH+3*H 
JW4=JW3+N 
Ft JW3,I4l=CVl(I,J) 

Ft JW3,I6)=-DV2(I,J) 

F  t  JH4, 1  5) =DV2 1 1 » J ) 

700  CONTINUE 

DO  800  J=1,N 

JW5=J+NW+2*N+3*H 

JW6=JW5+N 

Ft JH5,Ill J=CVltI,J| 

Ft  JW5,U3)=-DV2tI,J) 

Ft JH6,I12)=DV2tI,J) 

800  CONTINUE 
900  CONTINUE 


o  o  o  o 
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J0=NW+M+1 
J10=J0*2*M 
J20=J10+N 
I 30=3*M+1 
I40=I30+N 
F(JO, 11=1.0 
F(J10,I30)=1.0 
F { J2G, 1 40)  =  1 • 0 
DO  950  1=1, NW2 
GX( I ,1)=G( I I 
950  CONTINUE 

C  SOLVE  FOR  INITIAL  AND  FINAL  VALUE'.  OF  W 
CALL  SIMEQ  (F ,GX ,NH2, 1 ,K0) 

I  F{ KO.EQ.l )  GO  TO  3000 
DO  1000  1=1, NW 
W( I ,1 )=GX( 1,1) 

1000  CONTINUE 
C  FINAL  INTEGRATION 

CALL  INTEG  (GG, W, H, M,N, MM,NW, AX ,BX,CX ,DX, L ) 

C  INTEGRATION  DESTROYS  W(I,1I,  RESTORE  M(I,1) 

DO  1020  1=1, NW 
W( I ,1 )=GX( 1,1) 

C  CHECK  IF  THE  INTEGRATION  IS  ACCURATE 

I W=I+NW 

1020  AXU)  =  1.0-wn,MHi/GXUW,n 
AX  SHOULD  ME  MUCH  LESS  THAN  1 

IF  THE  ABOVE  CHECK  IS  UNSUCCESSFUL,  THE  STEP  SIZE  H  NEEDS 
TO  BE  DECREASED 

DETERMINE  SI  RAINS  AND  STRESSES 
DO  1500  J=1,NN 
Y=T*( J-1.0) /(NN-1.0J-T/2.0 
DO  1400  I  =1,MM 
EX=W(  (  M+l  I  *  I  )+Y*W(  (M«-2J  ,1) 

EY=W((3*M+2),I) 

EXY=W(2,I  HWI  <3*M+N+1),  I) +Y*W C ( 3*M+N*2 )  ,1  ) 

IFIM.LT.3)  GO  TO  1210 
DO  1200  11=3, M 
II1=II*M 

EX=EX*Y** ( 1 1-1 J ♦« (I  1 1 ,  II 
EXY=EXY*(II-1.C)*Y**I II-2)*W( I I, I ) 

1200  CONTINUE 

1210  IFIN.LT.3I  GO  TO  1310 
DO  1300  J J=3»N 
J  J3=3*M+J  J 
J J4=3*H+N+J J 

EY=EY+(  J  J-l  .0)  *Y*<«(  J  J-2  )*W(JJ3,IJ 
EXY=EXY+Y**(JJ-1)*WIJJ« ,1) 

1300  CONTINUE 
1310  SX=2.0*B*EX+BB*£Y 
SY=2.0*B*EY+BB*EX 
SXY=D*EXY 

SHEAR=0SQRT((SX-SY)4*2/4.0*SXY**2) 

SIGHAXI I , J»=SX 
SIGHAYI I, J J*SY 
TAUXYl I , JI=SXY 
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TAUMAXC I,J)=SHEAR 
1400  CONTINUE 
1500  CONTINUE 
3000  STOP 
END 

SUBROUTINE  INTEG  C GG, W,H, M,N, MM,NW, A, B, C, D, MM1 ) 

IMPLICIT  REAL*8CA-H,0-Z) 

DIMENSION  GGCMN»MN2) »WCNW,MM) ,A (NW) ,B( NW) ,C INW) ,DCNW) 

MX=1 

00  2000  KX=2,MM 
DO  1000  MMX=1 f MM1 
DO  10  I =1 , NW 
C  f 1 1=0.0 
10  CONTINUE 
K=-l 

100  K=K+1 

DO  200  1  =  1,  M 
I1=I+M 
I 2=I1+M 

D( I )=H*CWCI1,MXJ+CC III) 

D(  II )=H*(W( I2»MX)+CCI2)) 

DCI2)=0.0 
DO  120  JJ=1,M 
J J1=J J+M 

D { 1 2 ) =D ( 1 2 ) +H*GG ( I,JJ)*(W(JJl,MX)+C(JJl)) 

120  CONTINUE 

DO  150  JJ=1,N 

JJ1=JJ+M 

JJ2*JJ1+N 

JJ3=JJ+3*M 

JJ5=JJ3+2*N 

D(I2)=DCI2)+H*GGC I, JJI ) *C WC JJ3, MX )+CC JJ3) )+H*GGC I , JJ2)*(WC JJ5,MX) 
1+C ( J J5) ) 

150  CONTINUE 
200  CONTINUE 

DO  300  1=1, N 

I1=I+M 

I3=I+3*M 

I4=I3+N 

I 5=I4+N 

I6=I5+N 

D(I3)=Ht(w( I4»MX) +C  C 14) ) 

DC I4)=H*CWC I5»MX)+CC 15) ) 

DC I5)=H*(WC I6,MX)+CCI6)) 

DC  I 6)*0.0 
DO  220  JJ=1,M 
J J1*J J+M 

DC I6I=D! 161 +H+GGC II , JJ)+C WC JJ1,MX)+CC JJI) ) 

220  CONTINUE 

00  250  JJ*1,N 
J J1=J J+M 
JJ2=J  Jl  +  N 
JJ3-JJ+3+H 
JJ5»JJ3+2«=N 
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0<I6)=0<  !6MH*GG(  1 1 » J Ji )*( W ( JJ3 »MX l+C ( J J3I  )+H*GG(  1 1 » JJ2  )* 
l(W(JJ5»MXi+C(JJ5l) 

250  CONTINUE 
300  CONTINUE 

IFIK.EO.il  GO  TO  400 
IFtK.EQ.2l  GO  TO  500 
IF(K.EQ.3)  GO  TO  600 
DC  350  1=1, NW 

A<n*o(ii 
Cm*DU)/2.0 
350  CONTINUE 
GO  TO  100 
400  00  450  1=1, NW 
B( 1 ) -DC  I ) 

C ( I )=0( 1 1/2.0 
450  CONTINUE 
GO  TC  ICC 
500  DO  550  1*1, NW 
ClIhDII) 

550  CONTINUE 
GO  TO  100 
600  DO  700  1=1, NW 

W(  I,l|s«II,l)*U(IK2.0*BII)*2.0«C(I  )+D(  II  1/6.0 
700  CONTINUE 
1000  CONTINUE 

00  1100  1=1, NW 
1100  WU,KX)  =  WU,il 
2000  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  SIMEQ  ( A, B  ,NRCW,NRHS  ,K0 1 
C  THIS  SUBROUTINE  SOLVES  SIMULTANEOUS  EOLATIONS  WITH  REAL 

C  COEFFICIENTS 

IMPLICIT  RE AL*8  (A-H,C~Z) 

100  FORMAT!/'  SINGULAR  MATRIX'/) 

DIMENSION  A (NROW,NROW 1,B (NROW,NRHS I ,C (NW2) » 1 1 CNW2 1 , JJ  t  NW2, 2 
DO  1  1=1 ,NRCW 

1  11(11=0 

DO  17  I=i,NPOW 
D*O.DC 

DO  6  J=1 ,NROW 
IF(II(J)-i)  2,6,2 

2  DC  5  K=1 ,NROW 

IF( 1 1 CKl-1 1  3,5,21 

3  IF(CABS(C>-DABS(A(J,KI)I  4,4,5 

4  IR=J 
IC=K 

0=A{ J-KI 

5  CONTINUE 

6  CONTINUE 

II( ICI=II(IC1*1 
IF(IR-IC)  7,10,7 

7  DC  8  J=1 , NROW 

0=A ( IR , J  I 
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AIIR,J)=A<IC,J) 

8  A( IC, J)=D 

DO  9  J=1,NR HS 
D=B(IR,J) 

B(IR,J)  =  BUC,J) 

9  B< 1C* J)=Q 

10  JJ(I,1I*IR 
JJ(I*2)=IC 

C  ( 1  )  =  A( IC»IC) 

I F(CI I J )  11,21,11 

11  A(IC»IC)=1*D0 
DO  12  J*1,NR0W 

12  AUC»J)=A(IC»J)/C(I) 

DO  13  J=1,NRH$ 

13  BUCt  J)=B(IC,  Jl/C  (I  ) 

DO  17  J=l,NRdW 

I F( J-IC I  14,17,14 

14  D*A  ( J.  I C ) 

Al J,lCI=O.DO 
DO  15  K=1 , NROW 

15  A(  J,K)  =  A( J,K)-A(IC,KI*0 
DO  16  K*1 »NRHS 

16  B(J,K)*B<  J,K|-BUC,K)*D 

17  CONTINUE 

DO  20  1=1, NROW 
J=NR0W«-1-I 

IF( J J C J,1 l-J i  J , 2  ) )  18,20,18 

18  I R=J J ( J , 1 ) 

I  C  =  J  J  {  J  ,2  ) 

DO  19  J=1,NR0H 
D=A(J,IR| 

A(J,IR)=A(J,JCI 

19  A(J,IC)*D 

20  CONTINUE 
K0=0 
RETURN 

21  WRITE(6,100» 

KO=i 

RETURN 

END 


ooooonoo 
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APPENDIX  C 


FORTRAN  PROGRAM  BASED  ON  IBM  SYSTEM  360/75 
FOR  THE  IDENTICAL  ADHEREND  CASE 


Notation: 


E»E1  = 

0,01  = 

H»Hl  = 

B,B1  = 

A  * 

PtV,Cl,C2  •- 
MM* 

NN= 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  DETERMINE  STRESSES  IN  THE 
ADHESIVE  LAYER  OF  A  LAP  JOINT  WITH  IDENTICAL  ADHERENDS 
AND  GIVEN  LOADING- 

INPUT  *  E,El,Dt01,B,Bl,T,Tlf P, Vv Ml , A ,MM,NN 

OUTPUT  *  SIGMAX,  SIGMAY,  TAUXY,  T A UMAX 

IMPLICIT  REAL*8(A-E,H,0-Y),  C0MPLEX*16(F, G»  Z) 

C0MPLEX*16  CDEXP,  DCMPLX,  CDSQRT 

DIMENSION  S IGMAX(  MM,NN] , S I G MAY I  MM ,NN) , TAUXY (MM, NN) , TAUMAX(MM,NN) 
DIMENSION  UU12 ( MM ) »UU13( MM) ,VVOl ( MM) , VV02 ( MM) »  VV03 (MM) »  VV04I MM) 
DIMENSION  UUO(MM)  ,UU01  (MM)  ,UU02  (MM)  ,UU03(  MM),UU1(NM)»UUH(MM) 
DIMENSION  VVKMM)  ,VV1 1  <  MM ) ,  VV12  C  MM) ,  W13(  MM  )f  VV14(  MM) 

DIMENSION  F(5,5),G(5,1) 

DIMENSION  ZX(4),ZY(2),Z(11) 

C2=  CUV*A-P*(Tl*T) 

BT=B*T*2-0*8l*Tl 

CT=B*T/3-0*2-0*Bl*Tl 

AA1=B1*T1**2*T/(2.0*BT) 

AA2***H*T  /  (2  -  0*BT) 

AA3=Bl*Tl**2*(4.0/3.0*Tl*T)m*(Bl*Tl«‘*2«TW/2.0)  ) 
CV=T*CT*AA3-2.0*B1*T1**2 


E,  Ej  respectively 

G,  Gj  respectively 

X,  Xj  respectively 

Q/2  ,  Q^/2  respectively 

Overlap  Length,  a 

P,  V,  Mj,  M2  respectively 

Number  of  data  points  along x-axis,  (MM-1)  Ax  =  a 
Number  of  data  points  along  y-axis,  (NN-1)  Ay  =  t 
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CCC  =  Bl*Tl**3*T*(2.0*B*T+Bl*Tn 
CC--6.0*H*8T/CCC 

CC0=3.G*H*(2.0*BT-H*T)/CCC+12.0*Bt*D/CCC 

CC2=6.0*H*T*B1*T1**2/CCC+D*T*T/12.0+CC/H 

CC0=Bl*Ti**2-CT*T/2.0*AA3 

CC1=D*( AA3+1 .0) /CCD 

Q2=-D/ { 6. 0*CCD) 

Ql=3. 0*Q2 
Q0=6.0*Q2 
Q3=(CT*T/2.0)  /CCD 
Z(5)=P/(2.0*BT) 

FF=-CC1 

LI- 1 CC2/2.0 )**2-CCQ 
DO  100  N*1 » 2 

ZX(N)=CDSQRT<-CC2/2.0+(-l)**N*CDSQRT(ZZ)> 

ZX(N+2I =-ZX  ( N ) 

ZY(N)=(-1)**N*CDSQRT(FF) 

100  CONTINUE 

DO  200  J  =  l,4 

Ftl,  J)*( D*T*T/6. 04-2. 0*AA2*B1*T1 4*2 )*ZX( J)  +2.0*8 1*T  1**2* ( AA1- 
1T1*T/3.0)*ZX( J)**3 
F<2, J)=F( 1, J)*CDEXP< ZX( Jl+A) 

F(3»J)=B1*T1**2*( (Tl*T/3.0-AAl)*ZX( J) ++2-AA2) 

F(4»J)*F{ 3#J)*C0EXP(ZX{J)*A) 

200  CONTINUE 
Gll.l )*-V 
G(2,1)=V 

G*'3t  l>=-(C2-V*Al+P*Tl/2.0)/2.0*<1.0+AA2*CC/CC0)*Bl*Tl**2*Z(  5) 
G(4tl)=(Cl+V*Al)/2«0-P*Tl/4.Q+( 1. 0+AA2+CC/CC0 )*B1 *T1**2 *Z ( 5 ) 
CALL  SIMEQ  (F,G,4,1,K0I 
DU  300  N=*l,4 
Z(N)=G(N,1) 

300  CONTINUE 

DO  400  J=1 t 2 

F( 1, Ji^CV*ZY( J)**2 

F(2» J  >=F( 1 r JI*CDEXP(ZY( J!*A> 

F  (  3 »  J  J  =  ( D*T  *(AA3+1.0)*ZY(J)+B1*T1**2*(  AA3*T-4.0/3«  0*T1 )  *ZY(  J) 
I**3)*C0EXP<  ZY<J)*A) 

F(4,J)=(4.0/3.0*T1-AA3*T)+ZY( J)**2*CDEXP( ZY(J)*A) 

F  t  5 ♦ Jl  =  ZY( J ) 

400  CONTINUE 
F(l, 31=0.0 

F(2t3)=6.0+CV*Q2*A/CCl+CT*T*A 

F(3,3)=D*T*(AA3+1.0)/CCl*(3.0*Q2*A*A-6.0*Q2/CCl+Q3i+A*A/2.0*D*T 
UB1*TI»*2*{  (AA3*T-4.0/3.0*Tii*6.0*Q2/CCl  +  T) 

F(4»3)=A*( ( 8.0*T1-6.0*AA3*T)/CC1*Q2-T) 
F{5,3I=MQ3-6.0*Q2/CC1)/CC1 
F { lf4}*2.0*CV*Ql/CCl+CT*T 
F!2,4I*F(1,4) 

F<3,4)  =  D*T*A>M  ( AA3+U0)*2.0*Q1/CC1  +  1.0) 

F(4f4)=2.0*Ql*<4.0/3.0*Tl-AA3*T)/CCl-T 

F  ( 5  »4 )  =  0. 0 

F(l,5)=0.0 

F(2  »5)a0.0 

F( 3,5)=D4T*(1.0+Q04( AA3+1.0J/CC1) 
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F{4»5I=0.0 
F(5,5i=Q0/CCl 
G(l,l)=-P 
G(  2  »1  )  =  P 
G ( 3 » 1 )  =  V 

G(4,l  )=(CH-V*Al-P*Tl/2.0)/{  B1*T  1**25 
G(5,l)=0.0 

CALL  SIMEQ  (F,G,5,1,K01 
DO  500  N=6,10 
Z(N)=G( (N-5  )  #1 ) 

500  CONTINUE 

DO  800  N=1 » NN 

Y  =  T*(  N-1.0)/(NN-1.0)-T/2.0 

DO  700  M=1»MM 

X={ M~1.0)*A/(MM-1.0) 

V01=Z(  6)*ZY(1)*CDEXP( ZY(1 )*X)+Z (7)*ZY{2)*CDEXP(  ZY( 2 )*X>+3.0* 
1Z( 8)*X*X*Q2/CCl+<2.  *Q1*Z(9)*X*Q0*Z(10)*Q3*Z(8)-6.*Q2*Z<  8)/CCl 
25/CC1 

V02=Z«6  )*ZY<  15**2*CDEXP(  ZY( 1)*X  )  +  Z(  7)*ZY<  2)**2*CDFXP(  Z  Y  C  2  )  *  X )  ♦ 
16.0*Q2/CC1*Z(8)*X+2.0*Q1/CC1*Z(9) 

V03=ZI6)*ZY(1  )**3*CDEXP(ZYU)*X)  +  ZC  75  *Z  Y(  2  )  **3*CDEXP  (  ZY(2)*XH 
16 . 0*Q2/CC1*  Z ( 8  5 

V04=Z(6)*ZY(1) **4*CDEXP ( ZY  < 1)*X)  +  Z(75*ZY( 2 ) **4*CDEXP ( ZY(2 )*X) 

VV01 ( M) =V01 

VV02 ( M) =V02 

VV03 ( M ) =V03 

VVQ4 { M) =V04 

Ul=AA3*V01+Z(8)*X*X/2«0+Z(95*X+Z( 10) 

U11=AA3*V02«-Z(8)*X+Z(9> 

UU12(M)=AA3*V03+Z(8) 

UU13(M)=AA3*V04 
UU1 ( M )  =U1 
UU1 1 ( M) =U11 
F  V1=CC*  Z ( 5 ) /CCO 
FV1 1=0.0 
FV12=0. 0 
FV1 3=0. 0 
FV14=0«0 
DO  600  K=1 ♦ 4 

FV1  =  FV1«-Z<K)*CDEXP{  ZX(K)*X) 

FV11=FV1H-Z(K)*ZX(K)*CDEXP(  ZX(K)*X) 

FV12=FV12*ZIK)*ZX(K)**2*CDEXP(ZX(K)FX) 

FV13=FV13+Z(K)*ZX(K)**3*CDEXP(ZX(K)*X) 

FV 1 4=F V 14+Z  ( K ) *ZX ( K )  **4*CDEXP  I ZX { K) *X ) 

600  CONTINUE 
V i=FVl 
Vll=FVl 1 
V12=FV12 
VV1(M}=V1 
VV1 1(H) =V1 1 
VV12( M) =V12 
VV13 I M ) =F VI 3 
VV14( M) =FV14 
U01=AA1*V12*AA2*V1*Z( 5} 

UU01(M)~U01 
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UU02 ( M i =AA1*VV13( M)+AA2*Vll 
UU03 ( M 1= AA1 *VV14 ( M) +AA2*V12 
SX=2j0*B*(U01+Y*1'11  )+H*Vl 
SY=H*(U01+Y*U11)  *-2. 0*B*V1 
SXY=0*{  VC1*Y*V11+U1) 

SHEAR=DSQRT ( ( SX*SY  > ♦*2/4.0+SXY**2 ) 

S I GMAX ( M »  N) =SX 
S IGMAY ( M»  N) =SY 
TAUXY ( M » N )  =  SXY 
TAUMAX(M,N)=SHEAR 
STOP 
END 
C 

SUBROUTINE  SIMEG  (F , G, NROW, NRHS eKO) 

C  THIS  SUBROUTINE  SOLVES  SIMULTANEOUS  EQUATIONS  WITH  COMPLEX 

C  COEFFICIENTS 

IMPLICIT  C0MPLEX*16  <A-H,0-Z) 

REAL*8CDABS 

DIMENSION  F<5,5) ,G(5,l) ,C<5), 1 1 < 5 ) , JJ ( 5 ,2  ) 

IOO  FORMAT ( /*  SINGULAR  MATRIX'/) 

DO  1  1  =  1*  NROW 

1  IlllhO 

DO  17  1=1, NROW 
D= ( 0. 00,0.00) 

DO  6  J=1 , NROW 
IF{IUJ)-1)  2,6,2 

2  DO  5  K=1 » NROW 

IF(  II(K)-1 )  3,5,21 

3  I  F ( CDABS ( D)-CDABS(F( J,K) ) )  4,4,5 

4  I  R  =  J 
I  C=K 

D=  F  (  J  ,  K  ) 

5  CONTINUE 

6  CONTINUE 

II(  10=11(10  +  1 
IF(IR-IC)  7,10,7 

7  DO  8  J=1 , NROW 
D=F(  IR,  J) 

F(  IR  *  J  )=f  (  IC,J) 

8  F ( IC, J)=D 

DO  9  J=l, NRHS 
D*G( IR, J) 

G( IR, J)=G( IC, J) 

9  Gt IC* J)=D 
10  JJ( I , 1 ) =1 R 

JJ< I , 2) =IC 
C( I )=F ( IC, IC) 

I F ( CDABS( C(  I ) ) )  11,21,11 
II.  F<  IC,IC)  =  (1.D0,0.D0) 

DO  12  J  -l  ,NROW 

12  F ( IC, J )  =  F(  IC, J)/C (I  ) 

DO  13  J=1 , NRHS 

13  Gl IC, J)=G( IC, J)/C(I ) 

DO  17  J=l,NRGW 

I F ( J-!C  )  14,17,14 
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14  D=F(J,IC) 

F(J,IC)=(0.D0»0«D0) 

DO  15  K*1 , NROW 

15  F( J,K)*F( J,K)-FIIC,K1*D 
DO  16  K=1,NRHS 

16  G(  J,KJ-G< IC,K)*D 

17  CONTINUE 

DO  20  1*1, NROW 
J=NR0W+1-I 

IF( JJ( J,1 ))  18,20,18 

18  I R=J  J ( J , 1 ) 

IC  =  J J  C J,2) 

DO  19  J*1,NR0W 
D=  F { J ,  I R  ) 

FCJ,IR1  —  F(J,IC> 

19  F(J,IC)=0 

20  CONTINUE 
K0=0 
RETURN 

21  WR I TE ( 6 , 100  ) 

K0*1 

RETURN 

END 
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FIG.  I  Stress  due  tc  Stretching  of  Adherends 

(a)  Joint  Prior  to  Stretching,  (b)  Joint  after 
Stretching 


32 


t  ' 


(a) 


(b) 


FIG.  3  Lap  Joint  Under  a  General  loading 

(a)  Forces  Acting  on  Adherends  , (b) 
Forces  at  Joint  Edges 
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A  =  2"  for  Specimen  I 
A  =  3"  for  Specimen  2 


Fig.  6  Dimensions  of  Test  Models 


Photoelastic  Test  Bench 


Fringe  Order  at  Center  of  Disk 


Fig.  9  Calibration  of  Araldite  6020  Epoxy  to  Determine 
Fringe  Value 
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Fig.  10  Theoretical  Distribution  of  T,xy  at  Mid- plane  for 
Different  Values  of  m 

( Ei/E  =20  ,t,/t  =  4,  Va=  .5,  V«0,  n  *  I ) 


I 

l > 

I 

I 

I 

I 

I 


1? 


K- 


I 

t 


100 


Fig.  II  Theoretical  Distribution  of  Txy  at  Interfacial 
Plane  for  Different  Values  of  m 
(E|/e=20,  t,/t=4,  t|/a=.5,  V=0,  n  =  I ) 


0.2 


0.4 


0.6 


0.8 


1.0 


Distance  Along  Overlap  Length,  x/o 

Fig.  12  Theoretical  Distribution  of  0*y  at  Mid-plane  for 
Different  Values  of  m 

(E./E  *20,  t,/t*4,  ti/a*  .5,  V*0  ,  n*  I ) 


Normal  Stress  Factor,  Cy/T.av 
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Fig.  13  Theoretical  Distribution  of  O’ y  at  Interfacial 
plane  for  Different  Values  of  m 
(Ei/E  *20,  ti/t  *4,  h/a *  .5 ,  V*0,n*0 
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Fig.  15  Theoretical  Dietribution  of  O’ %  at  Interfacial 
Plane  for  Different  Values  of  m 


( Ei/E * 20  , ti/t »  4,  t./a  *  .5,  V*0,  n*  I) 
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X,  Distance  Along  The  Overlap  Length,  in. 


Fig. 17  Comparison  of  Theory  and  Experiment (y«  0 plane  ) 
a  « 2in.,  V/P  «  0 


Fringe  Ord 
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X,  Distance  Along  The  Overlap  Length,  in. 


Fig.  18  Comparison  of  Theory  and  Experiment  (y*  t/2  plane) 
a8  2 in. ,  V/P  3  0 


Fringe  Order 


108 


X  ,  Distance  Along  The  Overlap  Length,  in. 

Fig.  19  Comparison  of  Theory  and  Experiment  (y»0  plane) 
a  =  3in.,  V/P=  0 
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Fig. 23  Comparison  of  Theory  and  E xperiment (y» 0  plane ) 
a  ■  3in. ,  V/P  ■  0.  6 
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Fig.  24  Comparison  of  Theory  and  Experiment  (y»t/2plane) 
as  3in ,  V/P  «  0.6 
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Fig.  26  Comparison  of  Experimental  and 
Theoretical  Fringe  Patterns 
(q=  2  in.,  V/P  =  0.5,  0.6) 
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Experiment  (V/P  =  0.7) 


4.5  3.5 

Theory  (V/P  =  0.7) 
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Theory  (V/P  =  0.8) 


Fig.  27  Comparison  of  Experimental  and 
Theoretical  Fringe  Patterns 
(a  =  2  in.,  V/P  =  0.7,  0.8) 
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Experiment  ( V/P  =  0 ) 


3  5  4.5  3.5 


Theory  ( V/P  =  0  ) 


5.5 


Experiment  (V/P  =0.4) 


Fig.  28  Comparison  of  Experimental  and  Theoretical 
Fringe  Patterns  ( a  =  3  in.,  V/P  =  0.0,  0.4) 
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Experiment  {V/P  =  0.5) 


Theory  (  V/P  =  0.5  ) 
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Theory  (  V/P  =  0.6  ) 


Fig.  29  Comparison  of  Experimental  and  Theoretical 
Fringe  Patterns  (  a  =  3  in.,  V/P  =  0.5,  0.6) 
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Experiment  (V/P  =0.7) 
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“Theory  (V/P  =  07) 


Experiment  (V/P  =  0.8) 


Fig.  30  Comparison  of  Experimental  and 

Theoretical  Fringe  Patterns  (a  =  3in.,  V/P 


Shear  Stress  Concentration  Factor,  (Xxy  )max  /Xav 


Fig. 31  Shear  Stress  Concentration  Factor  at  the  Mid -plane 
of  Adhesive  Layer  as  a  Function  of  t./t,E./E ,  and  Sf, 
for  Vf  »0 


Parameter  ti/t 


Fig. 32  Shear  Stress  Concentration  Factor  at  the  Inter  facial 
Plane  as  a  Function  of  t,/t ,  E./E,  and  Sf.for  Vf  =0 
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Fig.  33  Shear  Stress  Concentration  Factor  at  the  Mid -plane  of 
the  Adhesive  Layer  as  Function  of  t«/t ,  E«/E ,  and  Sf, 
for  Vf  *  2 


Shear  Stress  Concentration  Factor ,  (‘Lxyjmox/Xav 


Fig.  34  Shear  Stress  Concentration  Factor  at  the  Interfacial 
Plane  os  a  Function  of  t«/t,  E,/E ,  and  Sf,  forVf  =  2 
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Fig.  35  Maximum  Shearing  Stress  Factor  as  a  Fui 
t./t.E./E,  and  Si,  for  Vf«0 


Maximum  Shearing  Stress  Factor,  Xmax/Xov 
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Fig.  36  Maximum  Shearing  Stress  Factor  as  a  Function  of  t,/t  , 
E,/E  ,  and  Sf,  for  Vf  =  2 


max 


Fig.  37  Shear  Stress  Concentration  Factor  at  the 
Midplane  As  a  Function  of  Sf  and  Vf 


Shear  Stress  Concentration  Factor, 
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Concentration  Factor.  (CTy )max/T.av  •  a/t 
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Fig.  40  Cy  Concentration  Factor  at  the  Mid -plane  as  a  Function  of 
ti/t,  E'/E  ,  Nf ,  and  Vf 


ncentration  Factor,  {(Ty )  max  / tav  •  a /t 
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Parameter  t./t 

Fig.  41  O’ y  Concentration  Factor  at  the  Interfacial  plane  as  c 
Function  of  t./t,  E. /E  ,  Nf  ,  and  Vf 


Parameter  ti/t 


Fig. 42  (J x  Concentration  Factor  at  the  Mid -plane  as  a  Function  of 
ti/t ,  E./E,  Nf ,  and  Vf 


Fig.  43  CFx  Concentration  Factor  at  the  Interfacial  plane  as  a 
Function  of  ti/t,  E</E,  Nf .  and  Vf 


Parameter  ti/t 


Fig.  44  Maximum  Principal  Stress  Factor  as  a  Function  of  t./t ,  E./E 
Nf , and  Vf 


Concentration  Factor 
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Parameter  Nf 

Fig.  48  CTX  Concentration  Factor  at  the  Intert'aciai 
Piane  As  a  Function  of  N,  and 


Parameter  Nf 

ig.  49  Maximum  Principal  Stress  Factor  As 


Function  of  N,  and  V< 
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F'g-  50  ay  Concentration  Factor  at  the  Mid -plane 
As  a  Function  of  (N,)2  and  v. 
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Fig.  52  0*  Concentration  Factor  at  the  Mid-plane 
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As  a  Function  of  (Nf)  and  Vf 
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Fig.  54  Maximum  Principal  Stress  Factor  As  a 
Function  of  (N,  )2  and  V, 


Shear  Stress  Factor 
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Fig.  55  IKy  Distribution  at  Interfacial  Plane  (y=  t/2) 
As  a  Function  of  Vf 
(E,/E  =  20,  t|/t  =  4,  tl/a=.3,  m=2) 


Normal  Stress  Factor 
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Distance  Along  Overlap  Length,  x/a 

Fig.  56  Oy  Distribution  at  Interfacial  Plane  (y=t/2) 
As  a  Function  of  Vf 
(E,/e  =  20,  t,/t=4,  t , /a  =  . 3 ,  m  =  2) 
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57  CJX  Distribution  at  Interfacial  Plane  (y=t/2) 
As  a  Function  of  Vf 
(E|/e  =  20,  t,/t  =  4 ,  t, /a  =  .3 ,  m  =  2) 


Top  Interface  (y  =  t/2) 
Midplane  (y-O) 

Bottom  Interface  (y  =  -t/2) 


Fig.  59  Txy  for  Joint  Loaded  in  Pure  Moment 
^ei/E=2®’  f i /t  =  ^ *  *i/a=-3»  nr»  =  2) 
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Ratio  of  Adherend  Thicknesses,  t2/t, 

Fig.  62  (TXy)max  'n  Non- Identical  Adherend  Joint  for  V  =  0 


Ratio  of  Maximum  Stresses, 
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Ratio  of  Maximum  Stresses, 
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Fig.  64  {Ox)max 


in  Non  -  identical  Adherend  Joint  for  V  =  0 


Ratio  of  Adherend  Thicknesses,  t2/t. 

Fig.  65  (Txy)max  in  Non- Identical  Adherend  Joint  for 


=  m2  =  o 


Ratio  of 


Fig.  66  (C7y  )max  m  Non-lde 


Thicknesses,  t2/t, 


Adherend  Joint  for  M,=Mz=0 


Ratio  of  Maximum  Stresses, 
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M,  =  M2=0 


Fig.  67  (CFx)max  in  Non- identical  Adherend  Joint  for 


